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ABSTRACT 

We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy 
distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample 
of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS 
photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the 
shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to 
selective absorption, Ry, can be determined with an uncertainty of about 0.1 for most stars in high-extinction 
regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar 
spatial distribution at low Galactic latitudes (\b\ < 30°). Our results are in good agreement with the extinction 
normalization given by the lSchlegel et al.l (119981 SFD) dust maps at high northern Galactic latitudes, but indi- 
cate that the SFD extinction map a ppears to be consis tently overestimated by about 20% in the southern sky, 
in agreement with recent study by Schl afly et alJ (2010). The constraints on the shape of the dust extinction 
curve across t he SDSS and 2MA SS bandpasses di s favor the reddening law of O'Donnell Q994), but support 
the models by lFitzpatrickl (Il999h and lCardelli et al.l (Il989h . For the latter, we find a ratio of the total to selective 
absorption to be TJy = 3 .0 ± 0. 1 (random)±0. 1 (systematic) over most of the high-latitude sky. At low Galactic 
latitudes (\b\ < 5°), we demonstrate that the SFD map cannot be reliably used to correct for extinction because 
most stars are embedded in dust, rather than behind it, as is the case at high Galactic latitudes. We analyze 
three-dimensional maps of the best-fit Ry and find that Ry = 3.1 cannot be ruled out in any of the ten SEGUE 
stripes at a precision level of ~ 0.1-0.2. Our best estimate for the intrinsic scatter of Ry in the regions probed 
by SEGUE stripes is ^ 0.2. We introduce a method for efficient selection of candidate red giant stars in the 
disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along 
the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly 
available in a user-friendly format. These data can be used for studies of stellar number density distribution, 
the distribution of dust properties, for selecting sources whose SED differs from SEDs for high-latitude main 
sequence stars, and for estimating distances to dust clouds and, in turn, to molecular gas clouds. 
Subject headings: methods: data analysis — stars: statistics — Galaxy: disk, stellar content, structure, inter- 
stellar medium 
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1. INTRODUCTION 

From our vantage point inside the disk of the Milky Way, 
we have a unique opportunity to study an L* spiral galaxy 
in great detail. By measuring and analyzing the properties of 
large numbers of individual stars, we can map the Milky Way 
in a nine-dimensional space spanned by the three spatial co- 
ordinates, three velocity components, and the three main stel- 
lar parameters - luminosity, effective temperature, and metal- 
licity. In a series of related studi es, we used data obtained 
by the Sloan Digital Sky Survey (I York et al .1120001) to study 
in detail the distribution of tens of mi l Uons of stars in this 
multi-dimensional space. In lJuricet al.1 (l2008l hereafter JOS) 
we examined the spatial distribution of stars in the Galaxy; 
in Ivezic et al. (2008a, hereafter 108) we extended our analy- 
sis to include the metallicity distribution; and in iBond et al] 
(120101 hereafter BIO) we investigated the distribution of stel- 
lar velocities. In Juric et al. (in prep) we estimate stellar lu- 
minosity functions for disk and halo stars, and describe an 
empirical Galaxy model and corresponding publicly available 
modelling code that encapsulate these SDSS-based results. 

All of the above studies were based on SDSS data at high 
Galactic latitudes {\b\ > 30°). Meanwhile, the second phase 
of SDSS has delivered imaging data for ten ^ 2.5° degree 
wide stripes (in SDSS terminology, two independent observ- 
ing runs produce two interleaving strips, which form a stripe, 
see Stoughton et al. 2002) that cross the Galactic plane (the 
so-called SEGUE data, see lYanny et al.il2009l) . At least in 
principle, these data can be used to extend the above analysis 
much closer to the mid-plane of the Galaxy, and to search for 
evidence of effects such as disk warp and disk flare. 

However, at low Galactic latitudes sampled by SEGUE 
data, there are severe problems with the interstellar dust 
extinction corrections. High-latitude SDSS data are typi- 
caUy corrected for i nterstellar extinction using maps from 
ISchlegel et all ([19981 hereafter SFD). When the full SFD ex- 
tinction correction is applied to low-latitude data, the result- 
ing color-magnitude and color-color diagrams have dramati- 
cally different morphology than those observed at high Galac- 
tic latitudes. Models developed by J08 suggest that these 
problems are predominantly due to the fact that stars are em- 
bedded in the dust layer, rather than behind it (the latter is 
an excellent approximation for most stars at high Galactic 
latitudes), and thus the SFD extinction value is an overesti- 
mate for most stars. This conclusion is also supported by 
other Galaxy mo dels, such as Besan gon (Robin et al . 2003!) 
and TRILEGAL (iGirardi et al.ll2005l) . Therefore, in order to 
fully exploit SEGUE data, both the intrinsic colors of a given 
star and the amount of dust extinction along the line-of-sight 
to the star have to be known. Distances to stars, which can be 
derived using appropriate photometric parallax relations (see 
108), would then enable mapping of the stellar spatial distri- 
bution. The interstellar medium (ISM) dust distribution and 
dust extin ction properties are interes t ing in their o wn right 
e.g., (e.g., iFitzpatrick & Massall2009l; lDrainell201 ll and ref- 
erences therein). An additional strong motivation for quanti- 
fying stellar and dust distribution close to the Galactic plane 
is to inform the planning of the Large Synoptic Survey Tele- 
scope (LSST) survey, which is c onsidering deep multi-band 
coverage of the Galactic planed ('I vezic et al.ll2008bi) . 

The amount of dust can be constrained by measuring dust 

See also Chapters 6 and 7 in the LSST Science Book available from 
www.lsst.org/lsst/scibook. 



extinction and/or reddening, typically at UV, optical and near- 
IR wavelengths, by measuring dust emission at far-IR wave- 
lengths, and by employing a tracer of interstellar medium 
(ISM), such as HI gas. For example, in their pioneering stud- 
ies in the late 1960s, Shane & Wirtanen used galaxy counts, 
and [Knapp & Kerr (1974) exploited a correlation between 
dust and HI column densities to infer the amount of dust ex- 
tinction. The most widely used contemporary dust map (SFD) 
is derived from observations of dust emission at 100 /im and 
240 /im, and has an angular resolution of ~6 arcmin (the tem- 
perature correction applied to IRAS 100 /im data is based on 
DIRBE 100 /im and 240 fim data, and has a lower angular 
resolution of ~ 1°; see SFD for more details). It has been 
found that the SFD map sometimes overestimates the dust col- 
umn by 20-30% when the dust extinction in the SDSS r band, 
A,. ^ 0.85Av, exceeds 0.5 mag (e.g., Arce & Goodman 1999). 
Such an error may be due to confusion of the background 
emission and that from point sources. A generic shortcoming 
of the far-IR emission-based methods is that they cannot pro- 
vide constraints on the three-dimensional distribution of dust; 
instead, only the total amount of dust along the line of sight 
to infinity is measured. In addition, the far-IR data provide 
no constraints for the wavelength dependence of extinction at 
UV, optical and near-IR wavelengths. 

With the availability of wide-angle digital sky surveys at 
optical and near-IR wavelengths, such as SDSS and 2MASS 
(see ^for more details), it is now possible to study the effects 
of dust exti nction using many tens of millions of sources. For 
example, ( Schlafly et al. 2010, hereafter Sch2010) utilized 
colors of blue stars, and Peek_&Grav_es (2010) utilized colors 
of passive red galaxies, to estimate errors in the SFD map at 
high Galactic latitudes. In both studies, the dust reddening 
is assumed constant within small sky patches, and the color 
distribution for a large number of sources from a given patch 
is used to infer the mean reddening (Peek & Graves dub this 
approach "standard crayon" method). Traditional dust red- 
dening estimation methods where the "true" color of a star 
is determined using spectroscopy we re extended to the ex- 
tensive SDSS spectroscopic dataset bv'Schla flv & Finkbeineii 
(2010) and Jones et al. (2011); they obtained results con- 
sistent with the above "standard crayon" methods. Studies 
of dust extinction with SDSS data are limited to Ay < 10; 
2MASS data alone can be used to trace dust up to Ay ^ 20 
using near-infrared color excess method (Lomba rdi & Alvell 
l200UlLombardi et al.l20Tll;lMaiewskrer al. 201 1), fliough es- 
timates of stellar distances are not as reliable as with SDSS 
data. 

In this work, we extend these studies to low Galactic lati- 
tudes where stars are embedded in dust, and also investigate 
whether optical and near-IR photometry are sufficient to con- 
strain the shape of the dust extinction curve. We estimate 
dust extinction along the line of sight to each detected star 
by simultaneously fitting its observed optical/IR spectral en- 
ergy distribution (SED) using an empirical library of intrin- 
sic reddening-free SEDs, a reddening curve described by the 
standard parameters: Ry =Ay/E{B-V), and the dust extinc- 
tion along the line of sight in the SDSS r band, A^. We first 
select a dust extinction model using high Galactic latitude 
data and another variation of the "standard crayon" method 
that incorporates the eight-band SDSS-2MASS photometry. 
Our SED fitting method that treats each star separately allows 
an estimation of the three-dimensional spatial distributions of 
both stars and dust. The dataset and methodology, including 
various tests of the adopted algorithm, are described in §2. 
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Results are analyzed in §3, and a preliminary investigation of 
the three-dimensional stellar count distribution and the distri- 
bution of dust properties is presented in §4. The main results 
are summarized and discussed in §5. 

2. DATA AND METHODOLOGY 

We first describe the data used in this work, and then dis- 
cuss methodology, including various tests of the adopted algo- 
rithm. All datasets used in this study are defined using SDSS 
imaging data for unresolved sources. Objects that are posi- 
tionally associated with 2MASS sources are a subset of the 
full SDSS sample. Although the SDSS-2MASS dataset is ex- 
pected to provide better performance than SDSS data alone 
when estimating dust properties and intrinsic stellar colors, 
we also consider the SDSS dataset alone (hereafter referred 
to as "only-SDSS") because it is effectively deeper (unless the 
dust extinction in the SDSS r band is larger than several mag- 
nitudes). We start by briefly describing the SDSS and 2MASS 
surveys. 

2.1. SDSS Survey 

The properties of the SDSS are documented in 
Fuku2ita et al. ( 1996); Gunn et al. ( 1998); H022 et al. (200j); 
Smith et all (I2002h: IStoughton et al.1 (l2002l). iPieretd] 
(I2003h; llyezic et all ( 120041) ; iTucker et al.1 (120061) and 
iGunn et al.l (12006b . In addition to its imaging survey data, 
SDSS has obtained well over half a minion stellar spectra, 
many as part of the Sloan Exte nsion for Galactic U nderstand- 
ing and Exploration (SEGUE: lYanny et al.l 12009). Here we 
only reiterate that the survey photometric catalogs are 95% 
complete to a depth of r 22, with photometry accurate to 
^0.02 mag (both absolute and rms error) for sources not 
limited by Poisson statistics. Sources with r < 20.5 have 
astrometric errors less than 0.1 arcsec per coordinate (rms; 
Ipier et al. 2003), and robust star/galaxy separation is achieved 
for r < 2 1 .5 (ILupton et al.li200l. 

The SDSS Data Release 7 (lAbazaiian et alll2009h . used in 
this work contains photometric and astrometric data for 357 
million unique object£3> detected in 11,663 sq. deg. About 
half of these objects are unresolved, and are dominated by 
stars (quasars contribute about 1%, see JOS). A full discussion 
of the p hotometric quality contr ol for the SEGUE scans is de- 
tailed in lAbazaiian et al.[(l2009h . Briefly, median reddening- 
free colors (Qgri and 2,-,v) were calculated for ea ch field using 
mag nitudes computed by both the SDSS photo ([Lupton et aLl 
I2OOII) and Pan-STARRS PS (Magn ier et alJl2010l) image pro- 
cessing pipelines, and the position and width of the locus of 
points (corresponding to the stellar main sequence) were com- 
puted. Fields within 15° of the Galactic plane had a wider dis- 
tribution ((TQ(photo,PS) 0.035, 0.027 mag) than fields out- 
side the plane (crQ(photo,PS) ^ 0.021, 0.20 mag). It can there- 
fore be inferred that (unsurprisingly) the photometric preci- 
sion in the plane is slightly poorer than that at higher lati- 
tudes. A more direct comparison is provided by magnitude 
differencing the photo and PS photometry. For stars with 
14 < u,g,r,i,z < 20, the median PSF magnitude difference 
was found to be 0.014 mag within the plane versus 0.010 mag 
outside the plane. 

2.2. 2MASS Survey 

The Two Micron All Sky Survey used two 1 .3 m telescopes 
to survey the entire sky in near-infrared light (ISkrutskie et alJ 

" For more details, see http://www.sdss.org/dr7/ 



Il997h . Each telescope had a camera with three 256 x 256 
arrays of HgCdTe detectors, and observed simultaneously in 
the / (1.25 /im), // (1.65 /im), and Ks (2.17 /im, hereafter K) 
bands. The detectors were sensitive to point sources brighter 
than about 1 mJy at the IQa level, corresponding to limiting 
magnitudes of 15.8, 15.1, and 14.3, respectively (Vega based; 
for corrections to AB magnitude scale see below). Point- 
source photometry is repeatable to better than 10% precision 
at these limiting magnitudes, and the astrometric uncertainty 
for these sources is less than 0.2". The 2MASS catalogs con- 
tain positional and photometric information for about 500 mil- 
lion point sources and 2 million extended sources. 

2.3. The Main-Sample Selection 

The main sample is selected from the SDSS Data Release 7 
using the following two main criteria: 

1. unique unresolved sources: objc_type=6, binary pro- 
cessing flags DEBLENDED_AS_MOVING, SATU- 
RATED, BLENDED, BRIGHT, and NODEBLEND 
must be false, parameter nCHILD=0, and 

2. the model r-band magnitudes (uncorrected for extinc- 
tion) must satisfy rMod < 21, 

These criteria yielded 73 million stars (for an SQL query used 
to select the main sample see Appendix A). The distribution 
of selected sources on the sky is shown in Figure [T] 

For isolated sources, the r < 2 1 condition ensures that pho- 
tometric errors are typically not larger than 0.05 mag (see 
Fig. 1 in Sesar et al., 2007). For sources with r < 19, the 
errors reach their systematic limit of ^0.02 mag. When re- 
ported errors are smaller than 0.02 mag, we reset them to 0.02 
mag to acc ount for expected photom etric zeropoint calibra- 
tion errors dPadmanabhan et al.ll2008l). T he behavior of best- 
fit x^tif distributions described in 33.1.1l iustifies the reset of 
errors. Errors can be much larger for sources in complex envi- 
ronments, and sometimes reported errors are unreliable (e.g., 
when sources are closer than 3 ", the photometri c errors are 
overestimated, see Figure 14 in lSesar et al.ll2008l) . If flie cat- 
aloged photometric error is larger than 0.5 mag in the griz 
bands, or larger than 1 .5 mag in the u band, that data point is 
not used in the analysis (formally, we reset the magnitudes to 
999.9 and their errors to 9999.9 in publicly available files, see 
Appendix B). 

2.4. SDSS-2MASS Subsample 

Following Co vey et al.l (I2007h . acceptable 2MASS sources 
must have 2MASS quality flags rd_flag == 222, bl_flag 
== 111, and cc^lag == 0, and selected 2MASS sources are 
positionally matched to SDSS sources with a distance cut- 
off of 1.5". The combined SDSS-2MASS catalog contains 
~23 million sources. The wavelength coverage of th e SDSS 
and 2M ASS bandpasses are shown in Fig. 3 in Fi nlator et alj 
(120001) ■ The distributions of SDSS-2MASS sources in vari- 
ous color-c olor and color -magnitude diagrams are discu ssed 
in detail by Finlator et all (2000) and lCovev et"aLl (l2007h . We 
emphasize that practically all sources in an SDSS-2MASS 
point source sample defined by a /T-band flux limit are suf- 
ficiently bright to be detected in aU other SDSS and 2MASS 
bands. For orientation, main sequence stars selected by the 
condition K < 14.3 are closer than approximately l-2kpc. 

Similarly to the treatment of SDSS photometry, for stars 
with reported errors in the J, H, and K bands greater than 
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0.5 mag, we reset magnitudes and errors to 999.9 and 9999.9, 
respectively. We also reset photometric errors to 0.02 mag 
when reported errors are below this hmit (systematic error s 
in 2MASS photometry are 0.02 mag; ISkrutslde eta^ll997^ . 
The Vega-based 2MASS photometry is translated to SDSS- 
like AB system following Finla tor et aL (2000.) . 

Jab = J2mass + 0-^9 (1) 

Hab = H2MASS +1-37 

Kab = K2MASS + 1 .84 

Note that these corrections have no impact on fitting and re- 
sults (because the same corrections are applied to models and 
observations and thus cancel, see below), but are convenient 
when visualizing SEDs. 

2.5. Model Assumptions and Fitting Procedures 

There are two empirical results that form the basis of our 
method. First, the stellar locus in the multi-dimensional 
color space spanned by SDSS and 2MASS colors is nearly 
one dimensional (because for most stars the effective tem- 
perature has much more effect on colors than other physi- 
cal parameters, such as age and metallicity). The locus po- 
sition reflects basic stellar physics and is so well defined 
that it has been used to test the quality of SDSS photometry 
(jivezic et al. 2004), as well as to calibrate new photometric 
data ( High et al.,.2009,) . 

Second, the shape of the dust extinction curve can be 
described as a one-parameter family, usually parametrized 
by gy^ Ay/EiB - V) (Cardelli et al. 1989; O'Donnell 1994^, 
iFiTzpatricld [l99l IMtzpatrick & Massal 120091) . Using this 
parametrization, extinction in an arbitrary photometric band- 
pass A is equal to 

Ax = Cx(Rv)A,-, (2) 

where Ar is extinction in the SDSS r band, and C\(Rv) de- 
scribes the shape of the extinction curv^^ Hence, the ob- 
served colors can be fit using only three free parameters: the 
position along the locus, Ry, and (eq.|2]is not the only way 
to "close" the system of equations; for a detailed discussion 
see Appendix B). Some caveats to this statement, such as the 
fact that not all unresolved sources are found along the locus 
(e.g., quasars and unresolved binary stars), and that even for 
fixed dust properties Ar and Ax depend on the source spectral 
energy distribution, are discussed in quantitative detail further 
below. We note that it is not mandatory to adopt an extinction 
curve parametrization given by eq. |2] For example, we could 
simply adopt the Ax values determined for high Galactic lati- 
tude regions by Schlafly et al. (2010). However, large dust ex- 
tinction observed at low Galactic latitudes offers a possibility 
to constrain the shape of the dust extinction curve, and eq. |2] 
provides a convenient one-pa rameter description that wo rks 
well in practice (but see also iFitzpatrick & Mas"sal 120091 for 
a different functional parametrization with two free parame- 
ters). 

A s imilar method was recently proposed by iBailer-JonesI 
(|201 1), where a strong prior is obtained from measured 
(trigonometric) distances and a requirement that stars must be 
consistent with stellar evolutionary track in the Hertzsprung- 
Russell diagram (as opposed to our constraint that stellar col- 
ors must be consistent with the stellar color locus). Such a 

The often used parametrization of dust extinction curve, k(X -V) = 
E(X-V)/E{B-V), is relat ed to Ca via k(\ -V) = Rv {CxAy/Av - 1); 
ISchlaflv & Finkbeineil j20TOl) give Av/Ar = l .200. 



prior has the advantage of being able to easily distinguish gi- 
ant stars from main sequence stars. Unfortunately, trigono- 
metric distances are not available for the vast majority of stars 
in our sample. 

2.5.1. Fitting Details 

The bes t-fit empirical stellar model from a library de- 
scribed in ^2.61 and the dust extinction according to a Cx{Rv) 
parametrization described in ^2.71 are found by minimizing 
X^df defined as 

where c"''' are observed adjacent (e.g., u-g, g-r, etc.) 
colors (A^ = 4 for only-SDSS dataset, and = 7 for SDSS- 
2MASS dataset). The number of fitted parameters is k = 3 
for all parameters, and k = 2 when a fixed value Ry = 3.1 is 
assumed (see below). 

The model colors are constructed using extinction- 
corrected magnitudes 

mr = mt-Ax, (4) 
with A = (ugriz[JHK]), resulting in 

c"'"' = c''\t)+[Cx2{Rv)-CM{Rv)] Ar. (5) 

Here Al and A2 correspond to two adjacent bandpasses which 
define colors c'"°"' and c''''. Hence, by minimizing x^^y, we 
obtain the best-fit values for three free parameters: /?v. A, , and 
the model library index, t (intrinsic stellar color, or position 
along the locus). Once these parameters are determined, the 
overall flux normalization (i.e. apparent magnitude offset) is 
determined by minimizing Xpj^ for the fixed best-fit model. 

We minimize x^^y by a brute force method. All 228 library 
SEDs (see ^2.6b are tried, with dust extinction values in the 
range < A, < 10 with 0.02 mag wide steps. This is not a 
very efficient method, but the runtime on a multi-processor 
machine was nevertheless much shorter, in both human and 
machine time, than post-fitting analysis of the results. 

We investigate the impact of Ry by producing two sets of 
best-fit t and A, . First, we fixed Ry = 3.1, and then allow Ry 
to vary in the range 1-8, with 0.1 wide steps. The results for 
the two cases are compared and analyzed in the next section. 

The errors, cr,, are computed from photometric errors 
quoted in catalogs, with a floor of 0.02 mag added in quadra- 
ture to account for plausible systematic errors (such as cal- 
ibration uncertainties), as well as for the finite locus width. 
In principle, c, could be varied with the trial library SED to 
account for the varying width of the stellar locus. We have 
not implemented this feature because it does not dominate the 
systematic errors. 

For a given Ry value (whether constant, or a grid value in 
the free Ry case), once the minimum y^, Xmin' located, 
an ellipse is fit to the section of the surface defined by 
X^ < Xmin'^^-^'^ (i-S-' within la deviation for 2 degrees of 
freedom): 

xHt,Ar\Rv) = a(t-t*f + b{t-t*)(Ar-A;) + c(Ar-A;)^ (6) 

were f is the model index, and t* and A* are the best-fit values 
corresponding to Xmin- Using the best-fit parameters a, b and 
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c, the (marginalized) model and A, errors can be computed 
from 



(7) 



a A 
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Note that the b coefficient controls the covariance between t* 
and A*. The surface around the best-fit t/A^ combination is 
described well by an ellipsoid, although this error ellipse ap- 
proximation breaks down far from the best-fit. The surface 
for stars with Xm/n > 200 is not fit with an ellipse and such 
stars, contributing less than two percent of the entire sample, 
are instead marked as bad fits. 

2.6. The Covey et al. Stellar SED Library 

ICovev et"an (l2007h quantified the main stellar locus in the 
ugrizJHK photometric system using a sample of ^600,000 
point sources detected by SDSS and 2MASS. They tabulated 
the locus position and width as a function of the g - i color, 
for 228 g-i values in the range -0.25 < < 4.50. We 
adopt this locus parametrization as our empirical SED library. 
Strictly speaking, this is not an exhaustive SED library that 
includes all possible combinations of effective temperature, 
metallicity and gravity, but rather a parameterization of the 
mean locus and its width in the multi-dimens ional SDSS- 
2MASS color space. We note that ICovey et all fcOO ?) used 
the SED map to correct SDSS and 2MASS photometry for 
interstellar dust extinction. Because they only studied high 
galactic latitude regions where typically A,. < 0.1, errors in 
derived locus parametrization due to errors in the SED map 
are at most 0.01 mag, and thus smaller or at most comparable 
to photometric calibration errors. 

This g - i parametrization reflects the fact that the stel- 
lar effective temperature, which by and large controls the 
g-i color, is more important than other physical parameters, 
such as age (gravity) and metallicity, in determining the over- 
all SED shape (for a related discussion and principal com- 
ponen t analysis of SDSS stellar spectra see iMcGurk et al.l 
120101) . The adopted g - i range includes the overwhelming 
majority of all unresolved SDSS sources, and approximately 
corresponds to MK spectral types from early A to late M. 
Due to 2MASS flux limits, the stellar sample analyzed by 
ICovey et al. (2007) does not include f aint blue stars (those 
with r>l6 for r < 0.6; see F'ig. 4 i n lFinlator e t al. 2000). 
Consequently, the ICovev et"an (l2007h locus corresponds to 
predominantly metal-rich main sequence stars ([Fe/H] > -1) 
because low-metallicity halo stars detected by SDSS are pre- 
dominantly faint and blue (see Eig. 3 in 108). According to 
Galfast model (J08), stars detected in SEGUE stripes are dom- 
inated by metal-rich main sequence stars, although we note 
that the fraction of red giant stars in SEGUE is expected to be 
much larger than observed by SDSS at high Galactic latitudes 
(- 20% vs. ~ 5%). 

The adopted model library cannot provide a good fit 
for SEDs of unreso lved pairs of white and red dwarfs 
jSmolcic et al.l I2004J). hot white dwarf s (Eisenstei n et al.l 
I2006h . and quasars (IRichards et al.ll200l1) . whose SEDs can 
differ from the adopted library by many tenths of a magni- 
tude. Systematic photometric discrepancies at the level of a 
few hundredths of a magnitude are also expected for K and 



M giants, especially in the u band jHelmi et al.ll2003l) . Sim- 
ilar M band discrepancies are expected for metal-poor main 
sequence stars (108). Nevertheless, all these populations to- 
gether never c ontribute more than ^ 20% of the full sample 
jFinlator et allllOOOt ljuridetar '2008). and m most cases can 
be recognized by their large values of xL/r l^^^t in prin- 
ciple, additional libraries appropriate for those other popu- 
lations can be used a posteriori to fit the observed SEDs of 
sources that have large Xmm when using SEDs of main se- 
quence stars. This additional analysis has not been attempted 
here, though our results represent the first necessary step: 
finding sources with large xlim- 

2.7. Parametrization of Dust Properties 

To implement the fitting method described in 32.51 the 
shape of the extinction curve {C\, see eq. |2]l must be charac- 
terized. C\ in the SDSS bands was initially computed (prior 
to the beginning of the survey, to enable spectroscopic tar- 
geting) using the standard parametr ization of the extinction 
curve tCardelh et al.lil989i: lo'DonneliilT994i) with Ry = 3.1. 
The resulting values (Ca=1.87, 1.38, 0.76, 0.54, with A = 
u,g,i,z) are commonly adopted to compute the extinction in 
the SDSS bands, together with Ar given by the SED map via 
Ar = 2.75E(B-V). 

A preliminary analysis of the position of the stellar locus 
in the SDSS-2MASS color space suggested that the above 
C\ values need to be slightly adjusted (M eyer et al.l 12005). 
Eurther support for this conclusion was recently presented by 
Sch2010. Here we revisit the Meyer et al. analysis using an 
improved SDSS photometric catalog from the so-called stripe 
82 regioi£3lvezic et al. (2007). SDSS photometry in this cat- 
alog is about twice as accurate as typical SDSS photometry 
due to averaging of many observations and various correc- 
tions for systematic errors. The SDSS-2MASS subset of that 
catalog includes 102,794 sources unresolved by SDSS (out of 
about a million in the full sample), which also have a 2MASS 
source with K < 14.3 within 1.5". The results of our analy- 
sis provide an updated set of C\ coefficients, which are then 
used to select a dust extinction model for generating the re- 
quired C\(Rv) dependen ce. Sim ilarly to a recent analysis by 
Sch2010, we find that the lO'Do nnell ( 1991) m odel can be re- 
jected , and adopt the CCM dust extinction law (ICardelli et al.l 
fl989h . 

2.7.1. Determination of the locus shifts 

The interstellar extinction reddens the stellar colors and 
shifts the position of the whole stellar locus at high Galac- 
tic latitudes, where practically all stars are located behind the 
dust screen. At high Galactic latitudes, distances to an over- 
whelming majority of stars are larger (> 100 pc) than the char- 
acteristic scale height of the interstellar dust layer (^^70 pc, 
J08). Both the amount of reddening and its wavelength de- 
pendence can be determined by measuring the locus position 
and comparing it to the locus position corresponding to a dust- 
free case. The latter can be determined in regions with very 
small extinction (Ar ~ 0.05) where errors in the SED extinc- 
tion map as large as 20% would still be negligible. 

We measure the locus position in the seven-dimensional 
SDSS-2MASS color space using an e xtended versi on of the 
"principal color" method developed bv llvezic et al.l (12004.) to 
track the quality of SDSS photometric calibration. We utilize 

^' Available from http://www.astro.washington.edu/users/ivezic/sdss/catalogs/stripe82.htn: 
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six independent two-dimensional projections spanned by the 
r-K and X-r colors, where X = u,g,i,z,J and H (see Fig.|2]i. 
Since the extinction in the 2MASS K band is small and fairly 
model and -independent (A^/A^ = 0.133 foiRy = 3.1, with 
only a ^ 1 0% variation over the range of plausible Ry and dust 
models, as discussed further below), the locus shifts in the 
r-K direction provide robust constraints for Ar- For exam- 
ple, a 10% uncertainty in the A^/A^ ratio results in only 1.5% 
uncertainty in Ar determined from a given Ar-Aj^ value. We 
determine these shifts iteratively, starting with Ar given by 
the SFD map, and adjusting Ar until the observed and cor- 
rected r-K color distributions agree in a maximum likeli- 
hood sense. This determination of A,- is very similar to the 
"blue tip" method introduced in Sch2010. The two main dif- 
ferences are due to the addition of 2MASS data. First, the 
low-metallicity faint blue stars are not included in the sam- 
ple analyzed here. Such stars could systematically influence 
the locus morphology and reddening estimates based on the 
"blue tip" method; nevertheless, our results are in good agree- 
ment with the Schlafly et al. results, as discussed below. Sec- 
ond, the availability of the K magnitudes enables a robust and 
straightforward determination of Ar, without any considera- 
tion of the SFD map. For a detailed discussion of these issues, 
see Appendix C. 

After Ar is estimated from the r-K color offsets, the lo- 
cus offsets in the A - r directions then provide constraints for 
the extinction wavelength dependence, Cx. We measure these 
offsets using principal colors. Pi and P2, with Pi parallel to 
the blue part of the stellar locus, and P2 perpendicular to it 
(see the top left panel in Fig.|2]for illustration of the principal 
axes, and for a comparison of the locus orientation with the 
direction of the standard reddening vectors). The blue part of 
the stellar locus at the probed faint magnitudes (14<r<17) 
includes mostly thick disk stars with distances of the order 1 
kpc or larger, which are thus beyond all the dust. 

We measure the position of the blue part of the locus in 
each A-r vs. r-ZT diagram using stars with 1.5 < r-K < 2.5 
(approximately; the range is enforced using the Pi (A) color). 
The blue part of the locus is parametrized as 



Pi ( A) = cos(0a) (r-K) + &m(9x) ( A - r) + ci( A) (9) 



and 



P2( A) = - siniOx) (r-K) + cos(0x) ( A - r) + C2(A). (10) 

The best-fit angle 9x found using stripe 82 dataset is equal 
to (61.85°, 33.07°, 14.57°, 23.47°, 34.04°, 43.35°). for 
A = (u,g,i,z,J,H). The values of ci and 02 are completely 
arbitrary; we set cx = 0, and determine C2(A) by requiring that 
the median value of P2(A) color is (c2=0.463, 0.434, 0.236, 
0.424, -0.048, -0.019, for u,g,i,z,J,H, respectively). Given 
the locus shift AP2(A), and Ar determined from the r-ZiT color 
offset (or alternatively from the APi offsets), the correspond- 
ing A A can be determined from 



Ca = — = l+tan(6lA)(l-— )+ — . 

Ar Ar COS(6x) Ar 



(11) 



Assuming a constant A;^ /Ar ratio, it is straightforward to com- 
pute the error of this estimate. 

The locus position must be measured over a sky area where 
the amount of dust and dust properties can be assumed con- 
stant. The smaller the area, the more robust is this assump- 
tion. However, the chosen area cannot be arbitrarily small 
because the error in the locus position, and thus the Cx error. 



is inversely proportional to the square root of the star counts. 
Within the analyzed stripe 82 region, the counts of SDSS- 
2MASS stars in the blue part of the stellar locus never drop be- 
low 70 stars deg~^. We bin the data using 4° wide bins of R. A. 
(with \Dec\ <1.27 deg., an area of ^10 deg^ per bin), which 
guarantees that random errors in Aa never exceed ^2% (even 
for the M band, and a factor of few smaller in other bands). 
In addition, we consider four larger regions: the high-latitude 
northern sky with b > 45°, split into Z < 180° and / > 180° 
subregions, a northern strip defined by 30° < b < 45°, and a 
southern strip defined by -45° < b < -30° (for these regions, 
we use SDSS DR7 photometry). 

2.7.2. Interpretation of the locus shifts and adopted dust extinction 

model 

We find that the variations in the shape of the extinction 
curve across the 28 R. A. bins from Stripe 82 region are consis- 
tent within measurement errors. The values of Cx obtained for 
the whole Stripe 82 region are listed in the first row in Table[T] 
Practically identical coefficients are obtained for the southern 
strip defined by -45° < b < -30°. The extinction curve val- 
ues for the northern sky are consistent with the southern sky, 
and we recommend the entries listed in the first row in Table\l\ 
for correcting SDSS and 2MASS photometry for interstellar 
dust extinction. One of the largest discrepancies is detected in 
a region from the northern strip defined by 30° <b <A5° and 
0° < / < 10°; and these values are listed in the second row 
in Table [1] Nevertheless, the north vs. south differences are 
not large, and, using models described below, correspond to 
an Ry variation of about 0. 1 . 

Much larger north vs. south differences are detected when 
comparing the best-fit Ar values to the SFD map values. The 
accuracy of the Ar determined here is about 3-10%, depend- 
ing on the amount of dust. We find that the SFD A,- values 
are consistently larger by about 20% than our values deter- 
mined across the southern hemisphere. Interestingly, no such 
discrepancy is detected across the northern sky, to within mea- 
surement errors of ~5%. In several isolated regions, the dis- 
crepancies are much larger. For example, in a region defined 
by -45° <b < -30° and 157° < / < 160°, the SFD values 
appear overestimated by 50% (the median value of A,- in that 
region given by the SFD map is 1.3). These results are simi- 
lar to those presented in Sch2010, where the spatial variation 
of errors in the SFD map and their possible causes are dis- 
cussed in more detail. The conclusion that the SFD A,- values 
are consistently overestimated in the southern hemisphere is 
als o consistent with the results based on galaxy count analysis 
by lYasuda et al.l (120071) . which is essentially an independent 
method. 

We adopt the Cx values determined for the Stripe 82 re- 
gion (the first row in Table [T]) to select a dust extinction law 
used in subsequent fitting of SEGUE data. Using the same 
assumptions and code as Sch2010, we compute dust extinc- 
tion curves for three popular models, and for three different 
input stellar spectral energy distributions. As can be seen in 
Figure[3] the differences between the models are much larger 
than the impact of different underlying spectra. 

A comparison of the observational constraints and model 
predictions is summarized in Figure ID Following Sch2010, 
we use ratios of the reddening values for this comparison. The 
differences in the extinction curve shape between the south- 
ern and northern sky determined here are similar to their dif- 
ferences from the Sch2010 results, and are consistent with 
estimated measurement uncertainties. The .O'Donnell (.19941) 



Milky Way Tomography IV: Dissecting Dust 



7 



model predicts unacceptable values of the {Ai -Ai)/(Ai-A,) 
ratio for all values of Ry- The other two models are in 
fair agreement with the data. Due to a slight offset of the 
Sch2010 measureme nts, they argued that the CCM model 
jCardelli et al.l ll989) is also unsatisfactory, although the dis- 
crepancy was not as large as in the case of the O'Donnell 
([094) reddening law. 

Although none of the models shows a perfect agreement 
with the data, the discrepancies are not large. To further illus- 
trate the constraints from different bands, we determine the 
best-fit Rv and its uncertainty in each band using the CCM 
model. If a model is acceptable, the constraints from dif- 
ferent bands must be statistically consistent. As shown in 
Figure |5] this is indeed the case, and we obtain the best-fit 
Rv = 3.01 ±0.05. The systematic error of this estimate, im- 
plied by the variation of the extinction curve shape across 
the analyzed regions, is about 0. 1 . The corresponding figure 
for the F99 ( Fitzpatrick ,1999) reddening l aw is similar, with 
the best-fit TJy = 3.30 ± 0.1, while for the lO'DonneUl (119941) 
model, Rv = 3.05 ±0.05. However, for the latter, the pre- 
dicted extinction in the / band is inconsistent with the rest of 
the bands at about la level (see Figure|6]l. This inconsistency 
is the main reason for rejecting the O'Donnell model both 
here and by Sch2010. A comparison between the CCM and 
O'Donnell laws is further illustrated in Figure |7j it appears 
that polynomial fitting adopted by both CCM and O'Donnell 
(to the 7''' and 8''' order, respectively) has caused wiggles 
whose integral over SDSS bandpasses is the largest in the / 
band. The predicted values of the extinction curve for all three 
models, using their individual best-fit values for Rv, are listed 
in Table n 

For the rest of our analysis, we generate C\(Rv) values us- 
ing the CCM law and an F star spectral energy distribution 
(7000 K). The adopted curves are shown in Figure |8] and a 
few representative values are listed in Table |2] For compar- 
ison, we also list C\ values suggested by Sch2010, and the 
values compute d using extinction cu rve parametrization pro- 
posed by Fitzp atrick & Massal (l2009h . 

2.8. Illustration of the Method and Fitting Degeneracies 

To summarize, we make two basic assumptions when ana- 
lyzing observed SEDs of low-latitude stars (SEGUE stripes). 
First, we assume that the median stell ar locus in SDSS a nd 
2MASS bandpasses, as quantified by Covev et aLl (120071) at 
high Galactic latitudes, is a good description of stellar col- 
ors at all Galactic latitudes. Second, we assume that the nor- 
malized dust extinction curve, Ax/Ar, can be described as a 
function of single parameter, Ry = Ay / E{B-V). Therefore, 
for a given set of measured colors, four in SDSS-only case, 
and seven in SDSS-2MASS case, we fit three free parame- 
ters: stellar model (position along the one-dimensional locus), 
f, dust amount. A,., and Ry- 

When the number of measured colors is small, when the 
color errors are large, or when the sampled wavelength range 
is not sufficiently wide, the best-fit solutions can be degener- 
ate. The main reason for this degeneracy is the similarity of 
the stellar locus orientation and the direction of the dust red- 
dening vector (see Figure |2]l. This degeneracy is especially 
strong for stars in the blue part of the locus {g-i < 1 ) and re- 
mains even when SDSS photometry is augmented by 2MASS 
photometry (a photometric band at a wavelength much shorter 
than the SDSS u band is needed to break this degeneracy). 

Figure |9] illustrates an example of degenerate solutions in 
the r-i vs. g-r color-color diagram, and how degeneracies 



are partially broken when the i-z color is added to the data. 
Because the direction of the reddening vector in the i-z vs. 
r-i color-color diagram is essentially independent of Ry, the 
measured r-i and i-z colors provide robust constraints for f 
and Ar, irrespective of Ry. The addition of the measured g - r 
color to r - / and i-z colors then constrains Ry. 

Since the stellar locus in the i-z and r-i color-color dia- 
gram and the reddening vector are not perpendicular, the co- 
variance between the best-fit t and A,- values does not van- 
ish. The addition of other bands, e.g. 2MASS bands to SDSS 
bands, alleviates this covariance, but not completely (and only 
slightly for blue stars). We quantify this effect using simulated 
observations, as described below. 

2.9. Tests of the Method 

To test the implementation of minimization algorithm, 
and to study the dependence of best-fit parameter uncertain- 
ties on photometric errors, the amount of extinction, and the 
intrinsic stellar color, we first perform relatively simple Monte 
Carlo simulations and analyze the resulting mock catalog 
based on realistic stellar and dust distributions, and photomet- 
ric error behavior. 

2.9. 1 . The Impact of Photometric Errors 

In the first test, we study the variation of best-fit parameters 
with photometric errors, where the latter are generated using 
Gaussian distribution and four different widths: 0.01, 0.02, 
0.04, and 0.08 mag. The dust extinction curve shape is fixed 
to = 3.1, and we only use SDSS photometry. The noiseless 
"observed" magnitudes for a fiducial star with intrinsic color 
^- / = 1 .95 (roughly at the "knee" of the stellar locus in the r- i 
vs. g-r color-color diagram) andAr= 1.5, are convolved with 
photometric noise generated independently for each band, and 
the resulting "observed" colors are used in fitting. The errors 
in best-fit models and A,, are illustrated in FiguresfTOlandfTTl 

The median errors in the best-fit stellar SED, parametrized 
by the g-i color, are about twice as large as the assumed pho- 
tometric errors. When photometric errors exceed about 0.05 
mag, the best-fit A^ distribution becomes bimodal, with the 
additional mode corresponding to a solution with a bluer star 
behind more dust. Therefore, even the addition of the red z 
passband is insuficient to break the stellar color-reddening de- 
generacy when the photometry is inaccurate (this conclusion 
remains true even when 2MASS bands are added). Our fitting 
results should thus be trusted only for stars sufficiently bright 
to have photometric errors smaller than about 0.05 mag in 
most bands. With this constraint, the formal best-fit errors are 
typically within 20% of the true errors. 

2.9.2. The Reddening vs. Intrinsic Stellar Color Degeneracy 

In the second test, we have investigated the covariance be- 
tween the best-fit model and A, values. Here again the dust 
extinction curve shape is fixed to Ry = 3.1. Figure [T2] shows 
the surfaces for a blue and a red star, and for two values 
of A,., when only SDSS bands are used in fitting and Gaussian 
noise with a = 0.02 mag is assumed for all bands. The best-fit 
model-Ar covariance is larger for the bluer star, in agreement 
with the behavior illustrated in Figure |9] (the angle between 
the reddening vector and the stellar locus is smaller for the 
blue part of the locus, than for the red part). The A^ vs. g-i 
covariance does not strongly depend on assumed A^. When 
the 2MASS bands are added, the morphology of the sur- 
face is essentially unchanged (recall that Ry was fixed in these 
tests). 
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These tests show that our implementation of the ^ mini- 
mization algorithm produces statistically correct results, and 
that the accuracy of SDSS and 2MASS photometry is suffi- 
cient (for most sources) to break degeneracy between the dust 
reddening and intrinsic stellar color in case of a fixed dust 
extinction curve (Ry = 3.1). Nevertheless, the best-fit results 
should be interpreted with caution when photometric errors 
exceed 0.05 mag, especially for intrinsically blue stars. 

2.9.3. Tests Based on a Realistic Galfast Mock Catalog 

To quantify the expected fidelity of our best-fit parameters, 
including Ry, for a realistic distribution of stellar colors, pho- 
tometric errors and dust extinction, we employ a mock cata- 
log produced by the Galfast code (Juric et al., in prep). Gal- 
fast is based on the Galactic structure model from JOS and in- 
cludes thin-disk, thick-disk, and halo components. The stel- 
lar populations considered here include main sequence and 
post-main sequence subgiant and giant stars. All other pop- 
ulations, such as blue horizontal branch stars, brown dwarfs, 
white dwarfs and quasars, are expected to contribute only a 
few percent of the total source count at low Galactic lati- 
tudes relevant here. SDSS and 2MASS photometry is gen- 
erated using the Covey et al. (2007) SED library (using the 
g-i color provided by Galfast). The photometri c errors are 
modeled using parametrization given by eq. 5 in Ivezic et aTl 
(|2008b), and the best-fit values for 5cr limiting depth derived 
using cataloged errors for SDSS and 2MASS data (for SDSS 
ugriz bands: 21.5, 23.0, 22.8, 22.6 and 20.5, respectively; 
for 2MASS 7///*: bands: 17.0, 16.0 and 15.5, on Vega scale). 
The dust extinction along the line-of-sight to each star is as- 
signed using the three-dim ensional dust distribution model of 
lAmdres & Lepind (120051) . The shape of the dust extinction 
curve is fixed to the CCM model values for Ry = 3.1. The 
normalization of the extinction for a given line of sight is de- 
termined by requiring a match to the SFD map at a fiducial 
distance of 100 kpc (that is, a complex dust distribution is re- 
tained in two out of three coordinates). 

The intrinsic absolute magnitude and color distribution of 
stars in the simulated low latitude {\b\ < 5°) sample is very 
different from distributions seen with high latitude samples. 
The two main differences are much bluer intrinsic color dis- 
tribution, and a much larger fraction of red giants in the low- 
latitude dataset. The origin of these differences is illustrated 
in the top two panels in Figure [T3] As shown in the top left 
panel, the simulated sample is dominated by stars with intrin- 
sic / < 1 .2, and includes a large fraction of red giants (40% 
with Mr < 2). These giants pass the r > 14 selection cut due 
to large dust extinction (Ar ^ 3 mag for giants in the simu- 
lated sample). At high Galactic latitudes, most red giants are 
brighter than SDSS saturation limit r ^ 14). 

The distributions of modeled stars in the color-magnitude 
and color-color magnitude diagrams closely match SDSS and 
2MASS data (for an illustration see Figure fl4l i. The much 
redder observed colors of stars in SEGUE stripes, compared 
to high-latitude sky, are reproduced with high fidelity. For ex- 
ample, the median g-i color for the SEGUE / 110° stripe 
moves from 1.0 at r ^ 16 to 1.7 at r ^ 21; only 2% of stars 
with r ^ 21 have g-i < 1. For comparison, at high Galac- 
tic latitudes, the median g - i color also becomes redder for 
fainter stars, but reaches a value of 1 at r ^ 19.5, or over 
3 mags fainter than at low Galactic latitudes. Although the 
two sets of diagrams are encouragingly similar, there a few 
detailed differences: the observed diagrams have more out- 
liers, and a few diagrams (e.g., J-K vs. i-z and i-z vs. 



r — i) imply different reddening vectors than used in simula- 
tions (Ry = 3.1). We discuss these differences in more detail 
in the next section. 

The resulting mock catalog is processed in exactly the same 
way as catalogs with observations. Note that the simulated 
photometry is generated with the same SED model library and 
dust extinction curve as used in fitting. We analyze four dif- 
ferent fitting methods: we use both only-SDSS (four colors) 
and SDSS-2MASS (seven colors) photometric data, and we 
consider both Ry = 3.1 (the true value) and Ry as a free fitted 
parameter Only stars with r < 20 and K < 13.9 (Vega) are 
used in analysis; this cutoff results in the median photometric 
errors of 0.02 mag in the r band and 0.04 mag in the K band 
(and 0.06 in the u band, which is the only band where errors 
exceed the K band errors). There are about 94,000 simulated 
stars that satisfy these criteria (the simulated area covers 25 
deg^). We first analyze the fitting results when Ry is fixed to 
its true value, and then extend our analysis to fitting results 
when Ry is a free parameter. 

When Ry is fixed, the obtained x^,df distributions closely 
resemble expected distributions for 2 and 5 degrees of free- 
dom, with slightly more objects in the tails. For example, 
86% and 93% of the sample are expected to have x^df < 2 
for only-SDSS and SDSS-2MASS cases, while we obtained 
73% and 80%. The latter fractions remain the same when the 
r band and K band limits are relaxed by 1 mag. For further 
analysis, we only use stars with xidf < ^■ 

The bottom left panel in Figure[r3]shows the distribution of 
simulated stars in the intrinsic apparent magnitude vs. color 
space, where we use only-SDSS best-fit intrinsic g-i color 
and correct "observed" r band magnitudes using the best-fit 
Ar. Its overall similarity with the top left panel is encouraging. 
The main difference is at the blue edge, g-i < 0.3, with about 
20% of stars having best-fit g-i color biased blue (simulated 
sample essentially does not include stars with g-i < 0.3 be- 
cause this is turn-off color for thick disk stars which con- 
tributes stars in that magnitude-color range). The same stars 
also have overestimated A,.. These biases are the result of the 
reddening-color degeneracy and could be mitigated by adopt- 
ing a strong prior such as removing SEDs with g-i < 0.3 from 
the SED model library. 

The bottom right panel compares the best-fit Ar to the in- 
put value. The best-fit Ar is systematically larger than the 
input values by about 10%. This overestimate is due to color- 
reddening degeneracy discussed above: when A,- is overes- 
timated, the best-fit stellar color is biased blue. When the 
full SDSS-2MASS dataset is used, the outliers seen in the 
bottom right panel in Figure [13] disappear, and the Ar bias 
is smaller by a few percent (10% vs. 15% for blue stars). 
Overall, there is no dramatic improvement resulting from the 
addition of 2MASS photometry (the root-mean-square scatter 
for the Ar difference decreases from 0.42 mag to 0.33 mag 
when 2MASS photometry is added). 

We find that the best-fit values based on only SDSS data are 
biased when the u band errors are large: A,- by 0.27 (true val- 
ues are smaller) and by 0.2 mag (bluer) for stars with u 
band errors of --0.1 mag. When the SDSS-2MASS dataset is 
used, both bias values fall to about 2/3 of only-SDSS values. 
Therefore, accurate u band photometry is crucial for obtain- 
ing accurate best-fit results. In order to minimize the effects 
of this bias, we further limit the sample to stars with u band er- 
rors below 0.05 mag. Unfortunately, only 40% of stars satisfy 
this cut. 



Milky Way Tomography IV: Dissecting Dust 



9 



The true errors in both stellar color and Ar (as determined 
by comparing the best-fit and true values) are about twice 
as large as marginalized errors computed using eqs. [T] and 
P both in case of only-SDSS and SDSS-2MASS fits. This 
increased scatter is probably due to color-A,. degeneracies 
and deviations of the maximum likelihood contours from a 
two-dimensional ellipse approximation: the errors in g-i 
color and Ar errors are strongly correlated with a slope of 
6{g-i)/d(Ar) ^ -0.65 (when this correlation is used to "cor- 
rect" the best-fit color by sliding them along this relation, 
the residuals are consistent with photometric errors; in other 
words, the entire "additional" color scatter is along this re- 
lation). The root-mean-square (rms) scatter for A, errors is 
0.42 mag and 0.33 mag for only-SDSS and SDSS-2MASS 
fits (20% and 16% for relative errors, i.e., errors normalized 
by true Ar), and the rms for g- i color errors are 0.29 mag and 
0.23 mag, respectively. We note that these errors are valid for 
individual stars, which suffer from the color-reddening degen- 
eracy. When the results are averaged in small pixels on the 
sky, the scatter is significantly smaller (because the spread of 
stars along the color-reddening degeneracy manifold is fairly 
symmetric). For example, the rms error for Ar in 0.2x0.2 
deg^ pixels decreases by a factor 3-4, to a level of about 5- 
10% (depending on the line of sight direction and the median 

Ar). 

2.9.4. "Free-Rv " Case 

The analysis of fits with Ry treated as a free parameter re- 
vealed that SDSS data alone are insufficient to reliably con- 
strain Ry, while SDSS-2MASS dataset produced good results. 
Figure [TS] compares the two resulting distributions of best-fit 
Rv (the input value is /Jy = 3.1). When SDSS-2MASS pho- 
tometry is used, Ry can be determined with a bias of < 0.1, 
and a precision (rms) of 0.10 when all stars from the simu- 
lated sample with Xp^// < 2 are considered. The Ry error is 
not correlated with stellar color, nor with distance; Ar is the 
only parameter that controls the Ry error. As expected, the Ry 
error increases for small A,.. A good practical limit is A^ > 1, 
which guarantees bias below 0.1 and an rms of at most 0.3. 
The Ry error decreases with Ar, and drops to 0.15 at A,- = 2 
and below 0.1 for A^ > 4. For A^ < 1, the precision of Ry esti- 
mate significantly deteriorates; for 0.5 < A,- < 0.7, the median 
best-fit Ry becomes biased to 3.2, with an rms of 0.5. 

Unsurprisingly, the Ry error is much larger when using only 
SDSS photometry; when considering all stars with x^^/ < 2, 
the best-fit Ry is biased to 3.3, with an rms of 1.2, rendering 
it practically useless. The main reason for this poor perfor- 
mance are the facts that three free parameters are constrained 
using only four colors, and that these three parameters are 
strongly degenerate. The SDSS-2MASS dataset shows su- 
perior performance when Ry is a free parameter not because 
2MASS data can constrain Ry (using CxiRy) parametrization 
employed here), but because 2MASS data better determine A,- 
and intrinsic stellar color, which gives more leverage to SDSS 
data (mostly the u and g band) to constrain Ry. 

As a result of this test, we conclude that only Ry estimates 
based on SDSS-2MASS dataset should be used, and those only 
for stars with x^p^f < 2 and A^ > 1 . 

2.9.5. "Dusty" Parallax Relation 

The analysis of the mock Galfast sample uncovered an in- 
teresting possibility for identifying candidate red giant stars 
in SEGUE stripes. Distinguishing red giant stars using only 



SDSS colors is difficult even at high Galactic latitudes (off- 
sets from the main sequence stellar locus are at most 0.02- 
0.03 mag; for more details see iHelmi et al.|[2003b . and seems 
futile at low Galactic latitudes. However, the best-fit Ar con- 
tains information about distance to a star, and this fact can be 
used for dwarf vs. giant star separation. 

After obtaining the best-fit intrinsic g-i color, we compute 
distance to each star using a photometric parallax relation ap- 
propriate for main sequence stars (108). For red giants, the 
resulting distances are grossly understimated (for example, a 
red giant star with g-i= 1 has Mr ~ 0, while main sequence 
stars with the same color have M, ~ 6, resulting in a distance 
ratio of ~ 15 for the same apparent magnitude). However, be- 
cause red giant stars are much more distant away than main 
sequence stars of the same color, their best-fit values of A,- are 
also on average significantly different. The latter difference 
is a consequence of the fact that A,, is proportional to the dust 
column along the line of sight, which in turn is roughly pro- 
portional to distance (although not exactly because the dust 
number density varies with position). 

These differences in the best-fit A,, vs. main-sequence dis- 
tance behavior between main sequence and red giant stars are 
illustrated in the top two panels in Figure[T6] The dashed lines 
mark the region in the Ar vs. distance diagram dominated by 
simulated stars with Mr < 3 (as illustrated in the bottom left 
panel). Red giant stars are found in the upper left corner of 
this diagram because their (main sequence) distances are too 
small given their A,-: it takes about 1 kpc of dust column to 
produce A r ^ 1 mag and thus stars with A ^ > 1 should be fur- 
ther than '-^l kpc. 

This separation of red giant and main sequence stars in the 
Ar VS. distance diagram can be elegantly summarized via a 
relation that we dub "dusty parallax". First, using the median 
best-fit Ar in narrow distance bins for stars with best-fit main 
sequence distances D < 0.5 kpc (see the blob discernible in 
the lower left corner), we obtained a linear relationship 

A,. = 1.06-^. (12) 
kpc 

The best-fit coefficient of 1.06 mag kpc"' is in good agree- 
ment with the coefficient corresponding to true A,- and dis- 
tance for stars with M,. > 5, 1.13 mag kpc"' (and implies that 
a similar algorithm can be applied to real data). This relation 
can be employed to estimate distance from the best-fit Ar for 
all stars, and in turn absolute magnitude Mr via "dusty paral- 
lax" relation 

Mf™ = r-5 logio(0.94Ar)-Ar- 10. (13) 

A comparison of true Mr and Mf^* is shown in the bottom 
right panel in Figure[T6] The root-mean-square scatter for the 
(M,--M™'^) difference is 1.2 mag. 

The coefficient from eq. [12] reflects the spatial dis- 
tribution of dust gener ated using a smooth model from 
lAmdres & Lepind (l2005h . In reality, localized clumps of dust 
will result in larger estimated distances and thus some main 
sequence stars will be misinterpreted as candidate red giants. 
Nevertheless, the precision of this relation seems sufficient to 
broadly separate red giant and main sequence stars using their 
best-fit g-i color and Ar. 

In many ways, this "dusty" parallax relation is similar to 
the reduced proper motion (RPM) method (f or a detailed dis- 
cussion, see Appendix B in lSesar et al.l2008h ; the main differ- 
ence is that RPM estimates distance using its relationship with 
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proper motion (assuming a fixed true tangential velocity, dis- 
tance is inversely proportional to proper motion), while DPR 
estimates distance using a relationship between dust extinc- 
tion and distance. We return to this relation and the selection 
of red giants when analyzing real data samples in the next 
section. 

To summarize this testing section, the analysis of simu- 
lated datasets has revealed important limitations of the best- 
fit results, mostly stemming from the finite photometric pre- 
cision of SDSS and 2MASS surveys. Most notably, the 
SDSS dataset alone does not have enough power to reliably 
constrain Ry, and only Ry estimates based on SDSS-2MASS 
dataset should be used, and those only for stars with x^^^j < 2 
and Ar>\. The tests based on a mock Galfast catalog also 
demonstrated that the fraction of red giant stars in low Galac- 
tic latitude samples is much larger than observed at high 
Galactic latitudes. These conclusions are important for the 
interpretation of results described in the next section. 

3. ANALYSIS OF THE RESULTS 

We apply the method described in the preceding Section 
(and summarized in 32.8b in four different ways. We fit sep- 
arately the full SDSS dataset (73 million sources) using only 
SDSS photometry, and the SDSS-2MASS subset (23 million 
sources) using both SDSS and 2MASS photometry. We first 
consider a fixed C\ extinction curve determined for Stripe 82 
region (the coefficients listed in the first row in Table [TJ, and 
refer to it hereafter as the "fixed Ry = 3.V' case (although the 
best-fit CCM model corresponds to Ry = 3.0 ±0.1). These 
fixed-/?!/ fits are obtained for the entire dataset, including high 
Galactic latitude regions where dust extinction is too small to 
reliably constrain the shape of the extinction curve (i.e., Ry) 
using data for individual stars. To investigate the variation 
of Ry in high-extinction and low Galactic latitude regions, we 
use the CCM C\ curves discussed in Section l2.7.2l (and shown 
in Figure Is). In this "free Ry" case, we only consider the ten 
SEGUE stripes limited to the latitude range \b\ < 30°, which 
include 37 million sources in the full SDSS dataset, and 10 
mil lion so urces in the SDSS-2MASS subset. As discussed 
in 32.9.41 the "free Ry'' results are only reliable when based 
on the full SDSS-2MASS photometric dataset. We include 
the "free Ry'' only-SDSS results in the public distribution for 
completeness, but do not discuss them further 

The resulting best-fit parameter set is rich in content and its 
full scientific exploitation is far beyond the scope of this pa- 
per The purpose of the preliminary analysis presented below 
is to illustrate the main results and to demonstrate their relia- 
bility, as well as to motivate further work by others - all the 
data and the best-fit parameters are made publicly available, 
as described in Appendix B. 

We first analyze "fixed Ry" fits, and compare results based 
on only-SDSS data with those obtained using the full SDSS- 
2MASS dataset. This comparison shows that both datasets 
result in similar best fits, which adequately explain the ob- 
served dust-reddened SEDs of most stars in the samples. The 
main conclusion derived from the "free Ry" fits is the lack 
of strong evidence for a significant overall departure from the 
canonical value of Ry = 3.1. 

3.1. Fixed Ry Case 

Two sets of results based on a fixed dust extinction curve 
("fixed Ry = 3.1" case) are compared: those based on the full 
SDSS-2MASS photometric dataset whose seven colors pro- 
vide better fitting constraints, and those for a larger and fainter 



only-SDSS sample which includes only four colors. We be- 
gin with a basic statistical analysis of the best-fit x^pdf distri- 
butions. 

3.1.1. The best-fit x^pd / distributions 

The distribution of the best-fit X^pdf' separately for low- 
extinction and high-extinction regions, and for low-SNR and 
high-SNR sources (bright and faint), is shown in Figure [17] 
As evident, there is no strong dependence of the shape of the 
best-fit x^df SNR. In low-extinction regions (top two pan- 
els) the obtained x^,df distributions closely resemble theoreti- 
cal xj,df distributions with 2 and 5 degrees of freedom. This 
agreement is not too surprising because the empirical model 
library was derived using the same dataset, and essentially 
demonstrates that cataloged photometric errors for SDSS and 
2MASS are reliable. 

In the high-extinction regions (although we discuss here 
only a single SEGUE stripe, we have verified that our con- 
clusions are valid for all ten stripes), the core of the observed 
X^,df distributions is still similar to theoretically expected dis- 
tributions (computed for Gaussian error distributions, and as- 
suming that SEDs of all stars in the sample are well de- 
scribed by the model library), but tails are more extended 
than in low-extinction high-latitude regions. For comparison, 
about 70% of a sample is expected to have x^df < 1 -2 (valid 
for the low number of degrees of freedom considered here), 
while we obtained about 50% for the observed distributions. 
The increased fraction of red giants at low Galactic latitudes, 
increased but unrecognized photometric errors (e.g., due to 
crowding), more complex dust extinction curve behavior than 
captured by the adopted CCM model, as well as increased 
metallicity of disk stars, may all contribute to the tails of the 
observed xj,df distributions. 

For further analysis, we use subsamples of stars with r < 19, 
K < 14 (Vega scale), and Xpdf < 2, unless noted otherwise. 
These criteria select stars with relatively small photometric 
errors (typically < 0.05 mag in most bands) and whose red- 
dened SEDs are well described by the model SED library and 
the CCM extinction curve. About 50-60% of stars in only- 
SDSS sample, and 70-80% stars in SDSS-2MASS subsample, 
are typically selected by the adopted xj,df < 2 cut (for theoret- 
ical x^df distributions with 2 and 5 degrees of freedom, 86% 
and 93% of stars would satisfy this x^,df cut). 

3.1.2. The Northern Galactic Cap Region 

Due to small Ar for the b > 30° sky region (the median A,- 
from the SFD map is ~0.08 mag), the errors for best-fit A,- for 
individual stars can be as large as best-fit A,- itself when using 
only-SDSS fits {fixed Ry case). Both the formal Ar errors, and 
the differences between best-fit and SFD values for A^ begin 
to increase rapidly for r > 18 and become unreliable for r > 
19. This behavior is in agreement with tests described in 32.91 
and the behavior of SDSS photometric errors as a function 
of magnitude (even for blue stars, the median u band error 
is already 0.05 mag at r = 19, and 0.2 mag when stars of all 
colors are considered). 

Nevertheless, by taking a median value for typically several 
hundred stars per '^l deg^, a map can be constructed that re- 
produces the features seen in the SFD map (see the top left 
panel in Figure [TsTl. Quantitative analysis of the median dif- 
ferences between the best-fit A^ and the SFD A^ values shows 
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that the former are larger by about 50% on average, with a 
scatter of about 20%. This bias is probably due to color- 
reddening degeneracy and small extinction at high Galactic 
latitudes which is only a factor 2-3 larger than photometric 
errors. An additional effect contributing to this bias are zero- 
point calibration errors in SDSS photometry: the median dif- 
ferences between the best-fit A,- and the SFD Ar values show 
a structure reminiscent of the SDSS scanning pattern (see the 
top right panel in Figure [TSll. These coherent residuals imply 
problems with the transfer of SDSS photometric zeropoints 
across the sky. 

The median differences between observed and best-fit 
model magnitudes show deviations of up to 0.01 mag, and 
are largest in the / band, as illustrated in the bottom left panel 
in Figure [TsTl. Therefore, these relatively small local calibra- 
tion errors (each of the six scanning strips in an SDSS scan, 
i.e., the "camera columns", is independently calibrated) are 
mis-interpreted as a local extinction variation at the level of a 
few times 0.01 mag. 

With the addition of 2MASS photometry, the agreement 
with the SFD map improves. The best-fit Ar values are over- 
estimated, relative to SFD values, by only ^0.02 mag (25% 
on average), and the median differences do not show structure 
resembling the SDSS scanning pattern (see the bottom right 
panel in Figure [TsTl. We note that r < 18 selection limit (and 
A' < 14 in 2MASS case) results in about one star per the res- 
olution element of SFD map. Therefore, to significantly im- 
prove the spatial resolution of extinction map at high Galac- 
tic latitudes, a sample several magnitudes deeper than SDSS- 
2MASS sample is required. 

3.1.3. The SEGUE Stripes 

The main goal of this work is to determine extinction at 
low Galactic latitudes. We consider ten ^ 2.5°. wide SEGUE 
stripes with \b\ < 30°. The full SDSS sample includes 37 
million sources, with 10 million sources in the SDSS-2MASS 
subset. We find that results based on the two datasets are sim- 
ilar, though the latter is expected to produce more reliable re- 
sults. We first illustrate the behavior of best-fit A ^ as a function 
of distance for all stripes, and then provide more quantitative 
discussion of the differences in best-fit results in the next sec- 
tion, which is focused on a single stripe (/ ^ 110°). We also 
provide a comparison to the SFD extinction maps further be- 
low. 

A visual summary of the best-fit A,- using only-SDSS fits 
for the ten SEGUE stripes, in the range \b\ < 5° and for three 
distance slices ranging from 0.3 kpc to 2.5 kpc, is shown in 
Figure[T9] Distances to stars are determined by assuming that 
all sources are main sequence stars, and using photometric 
parallax relation from 108 with [Fe/H]=-0.4 (with the best-fit 
intrinsic colors). An expected scatter in metallicity of 0.2-0.3 
dex for disk stars corresponds to about 10-15% uncertainty in 
distance. Although not all sources are main sequence stars 
(such as red giants, which have grossly underestimated dis- 
tances, see ^3. 1.5l below for discussion), the fraction of main 
sequence stars in the samples is sufficiently large that the me- 
dian Ar is not strongly biased. Furthermore, sources whose 
SEDs are significantly different from the main sequence SEDs 
are not included: the figures are constructed only with sources 
that have the best-fit x^pdf < 2. We also excluded red giant 
candidates, as described below. 

It is easily discernible from Figure [19] that the extinction 
along the Une-of-sight (that is, Ar) increases with distance. 



On average, the stripes towards the Galactic center have more 
large-extinction (A,. > 1) regions. In several directions. A,- ex- 
ceeds several magnitudes and practically no stars are detected 
by SDSS. 

3. 1 .4. Selection Function Dijferences for only-SDSS and 
SDSS-2MASS subsamples 

Another projection of the sky position-distance-A^ space is 
shown in Figures|20](only-SDSS case) andl2T](SDSS-2MASS 
case). As evident, the morphology of these A,, vs. distance di- 
agrams differs significantly between the two subsamples. The 
main reason for these changes is different sample selection 
functions in the flux-color space - and not differences in the 
best-fit Ar or distance values which agree well on a star by star 
basis (see the next section). 

For only-SDSS case, the main selection criterion (in addi- 
tion to Xpdf < 2 in both cases) is r < 19 and the u band error 
limit of 0.05 mag. The latter condition is necessary to assure 
reliable fitting results when only four colors are used, and re- 
sults in a strong bias towards the blue end of the observed 
color distribution. In SDSS-2MASS case, it is sufficient to re- 
quire K < 15 (Vega) to obtain reliable fitting results because 
there are seven colors, and because this condition limits the K 
band and u band errors to about 0. 1 mag (with much smaller 
errors in other bands). This selection condition results in a 
strong bias towards the red end of the observed color distribu- 
tion. Due to their selection functions, the effective r band lim- 
iting magnitude for reliable only-SDSS samples varies from 
r - 17 at / = 1 to r - 15 at g-i = 3, while for SDSS- 
2MASS samples it varies from r ^ 17 at g - / = 1 to r ^ 20 
at / = 3. As a result, SDSS-2MASS samples contain many 
more nearby red dwarfs at distances below 500 pc, while only- 
SDSS sample extends further than SDSS-2MASS sample, to 
about 2.5 kpc. On average, about twice as many stars survive 
the quaUty cuts for SDSS-2MASS sample as for only-SDSS 
sample (although the latter typically contains about four times 
as many stars at \b\ < 5° before any selection). 

In the Ar vs. distance diagram, the selection cutoff for 
SDSS-2MASS sample is nearly vertical, and limits the sam- 
ple to distances below about 1.5 kpc (assuming Ar < 5 and 
main sequence stars). For only-SDSS sample, the upper limit 
on u band error introduces a diagonal selection boundary that 
excludes stars in the upper right corner With the selection 
criteria adopted above, the sample becomes limited to Ar < 2 
at a distance of about 1 kpc, with an overall distance limit of 
about 2.5 kpc. 

The slopes of Ar vs. distance relations along main-sequence 
locus seen in Figures |20] and |2T| constrains the local (within 1 
kpc) extinction per unit distance normalization to the range 
Ar/D = 0.7- 1.4 mag kpc"', with larger values correspond- 
ing to smaller angular distances from the Galactic center The 
variation of this normalization with Galactic longitude is con- 
sistent with the exponential scale length for thin disk stars 
obtained by J08 (Li = 2.6 ± 0.5 kpc, see their Table 10). Nev- 
ertheless, the variation fo Ar/D with longitude observed here 
is more complex than predicted by simple axially symmetric 
dust distribution model. 

3.1.5. The Selection of Candidate Red Giant Stars 

The Ar vs. distance diagrams based on SDSS-2MASS data 
(see Figure l2ni show an excess of sources in the top left cor- 
ner (the effect is not as strong for only-SDSS case because the 
selection effects due to the u band error limit, discussed in the 
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previous section, remove mo st of th ese sources). Based on a 
mock catalog discussion in ^2.9.51 these sources are consis- 
tent with red giant stars. Informed by their distribution, and 
clear separation from the locus of main sequence stars, we 
adopted the following criteria for the selection of candidate 
red giants: 

1. Best-fit main sequence distance below 1 kpc, Dkpc < 1, 

2. Best-fit extinction. A, > 1.5+ LSDi^pc, and 

3. Best-fit intrinsic color, 0.4 < g-i < 1.4. 

The first two criteria are based on the morphology observed in 
the A, vs. distance diagrams, and the third criterion removes 
outliers whose best-fit intrinsic colors are inconsistent with 
the color distribution for the majority of sources selected by 
the first two criteria. 

We applied these criteria to all ten SEGUE stripes and 
found that the fraction of selected stars varies significantly 
with Galactic longitude, from '^15% for stripes at / = 50° 
and I = 70° to ^ 2% for stripes within 20°. from the Galactic 
anticenter The inclination of the main sequence stellar locus 
in the A^ vs. distance diagrams also varies with Galactic lon- 
gitude, with its slope (determined for distances up to 1 kpc) 
decreasing from about 2.0 mag kpc"' for the I = 50° stripe to 
0.6 mag kpc"' for the / = 187° stripe. Hence, our selection 
criterion #2 above could be improved by taking this variation 
into account (for the same reason, the proportionality "con- 
stant" in eq. [T2lvaries with longitude). 

The observed variation of the fraction of candidate red gi- 
ants with Galactic longitude represents a strong constraint for 
the Galactic structure models, and the change of A,, vs. dis- 
tance slope reflects the variation of dust number volume den- 
sity in the Galactic disk. Hence, the data presented here can 
be used to improve Galactic stellar population models such as 
Galfa st and TRILEGAL , and dust distribution models, such 
as the lAmores & Lepind (12005) model employed by Galfast. 
The required detailed analysis is beyond the scope of this 
work. 

3.1.6. Detailed Analysis of the I ~ 110° SEGUE stripe 

For a detailed analysis of the best-fit results, we select a 
single fiducial SEGUE stripe with Z '--^ 110°. A simple but 
far-reaching conclusion of the work presented here is that fits 
to intrinsic stellar SEDs and dust extinction on per star basis 
are capable of reproducing the morphology of observed color- 
color diagrams in highly dust-extincted regions. This success 
is illustrated in Figure]^ where six characteristic color-color 
diagrams constructed with observed SDSS-2MASS photome- 
try are contrasted with analogous diagrams constructed using 
best-fit results. We reiterate that the observed morphology 
in these diagrams at low Galactic latitudes is vastly different 
than at high latitudes (the latter is illustrated in the figure by 
the Covey et al. locus). 

When considering SDSS-2MASS sample, fits based on the 
full seven-color set and those restricted to the four SDSS col- 
ors produce quantitatively similar, though not identical re- 
sults. The root-mean-square (rms) scatter of the difference 
in best-fit intrinsic colors is 0.04 mag, and rms for best-fit Ar 
difference is 0.07 (the median A, is 1.9). For A,. ^ 5, the val- 
ues based on only-SDSS photometry become biased (larger) 
by about 4% relative to SDSS-2MASS values. A star-by-star 
comparison presented in Figure |23] shows a few regions (e.g. 



g-i ^ 1.5 and small g-i) where results can differ substan- 
tially; nevertheless, the fraction of affected sources is small 
and negligible when results are averaged over many stars. The 
latter point is illustrated in Figure |24l which compares the 
two Ar maps for stars at a limited range of distances. The 
two maps agree to better than 0.05 mag even in regions where 
Ar > 4. This agreement demonstrates that SDSS data alone 
are sufficient to obtain the best-fit intrinsic color and extinc- 
tion along the line-of-sight for the majority of stars (when Ry 
is fixed). In the rest of analysis we use SDSS-2MASS results, 
except in a few cases where we explore distances beyond 2 
kpc. 

A cross-section of the three-dimensional A,- map, based 
on only-SDSS sample from the / ^ 110° is shown in Fig- 
ure |25] As evident, the best-fit Ar increases with the stellar 
distance between 0.3 kpc and 2.5 kpc. It is noteworthy that 
the two quantities are determined independently (distance is 
computed a posteriori, from the best-fit apparent magnitude). 
A closer look at distances below 1 kpc using SDSS-2MASS 
dataset is shown in Figure |26] An impressive feature is the 
abrupt jump in Ar towards b ^ 2° for stars with distances 
above 0.9 kpc, thus providing a robust and fairly precise lower 
distance limit for that dust cloud! 

Differences between best-fit Ar values determined here and 
the SFD map are illustrated in Figure|27] Since the latter cor- 
responds to extinction along the line of sight to infinity, our 
values are systematically smaller in regions with large A, and 
similar at large Galactic latitudes, as expected. A detailed 
analysis of these Ar differences, when combined with stellar 
distance estimates, can provide valuable constraints for vari- 
ous ISM studies. For example, in Figure |28] we demonstrate 
good correspondence between the A, differences and the dis- 
tribution of molecular (CO) emission; our results imply that 
those molecular clouds must be more distant than 1 kpc, and 
that the substructure seen around h -2.5° is more distant 
that the one atb^ 2 ° (se e also Figure|26l and a more quanti- 
tative discussion in 34.11 ). Other SEGUE strips contain more 
examples where such "bracketing" of distances to molecular 
clouds can be attempted. 

We note that the SFD map is expected to sometimes fail 
at low Galactic latitudes not just because of stars being em- 
bedded in dust, but also because its construction relied upon 
accurate point source subtraction (which is only performed 
for 1^1 < 4.7°) and dust having a single temperature along 
each line of sight. These assumptions might be violated in 
the Galactic plane and may be responsible for regions where 
SFD values are smaller than in our maps (e.g., atb^ 10-20° 
in Figure 26). 

3.2. Free Ry Case 

If there is a significant discrepancy between the shape of 
assumed CCM extinction curve for Ry = 3.1 and that required 
by SEGUE data, photometric residuals between observed and 
best-fit magnitudes should show a correlation with best-fit 
Ar. Indeed, the failure to pass this test has revealed that our 
first instance of fitting erroneously used the O'Donnell extinc- 
tion curve (due to an error in "metadata management"). In 
this case, the photometric residuals (data "minus" model) in 
the / band showed a highly statistically significant correlation 
A/ = -0.015 Ar, which implied that the adopted C\ value in the 
/ band was too large by 0.015 (the results for other bands did 
not require a change of C\). The analysis of used C\ values 
clearly placed them on top of the O'Donnell model in the right 
panel in Figure|4] while the revised value moved the constraint 
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towards the CCM model curves. After our second fitting it- 
eration that correctly incorporated the CCM model, we re- 
gressed photometric residuals and best-fit Ar again and found 
much smaller residuals: A/ = -0.005 Ar and Az = 0.003 A,.. 
For no other bands were the slopes larger than statistical mea- 
surement errors of at most 0.001. These two relatively small 
corrections of C\ in the / and z bands result in a shift in the 
right panel in Figme^away from the CCM model curves, and 
to a point between the constraints obtained using stellar locus 
method for Stripe 82 and the northern Galactic hemisphere! 
That is, the required C\ modifications cannot be accomplished 
by adopting a CCM model curve for a different Ry (nor using 
any of the other two considered models). Hence, SEGUE data 
"knew" that (independent) empirical constraints on the shape 
of dust extinction curve from the high-latitude sky are better 
than the CCM model for = 3 . 1 ! 

The above analysis of photometric residuals shows that 
there is no a priori reason to expect a significant departure 
from the canonical Ry = 3.1 value when Ry is considered a 
free fitting parameter. Nevertheless, it is possible that local- 
ized regions in the Galactic disk have a different Ry distribu- 
tion, and given the unique nature of our sample, such a study 
is worthwhile. Th e analy sis of fitting results for a mock cat- 
alog described in ^Z93l showed that only the SDSS-2MASS 
dataset can be expected to provide useful constraints on Ry, 
and this is the fitting case analyzed here (for completeness, 
public data distribution includes also only-SDSS case). 

A comparison of best-fit intrinsic colors and Ar between 
fixed-/?!/ and free-/?;/ cases is shown in Figure |29] While 
for some sources results can differ substantially, the fraction 
of discrepant sources is small. The resulting distribution of 
sources in the Ar vs. distance diagram, shown in Figure |30] 
for fiee-Ry case, is similar to that based on fixed-Ry case 
(compare to Figure l2Tli. A comparison of best-fits results for 
fixed-/?!/ and free-Ry cases shown in Figure[3T|reveals that fit 
residuals are not significantly smaller when Ry is free. 

The median Ry, as a function of the position in the A,- vs. 
distance diagram, is shown in Figure [32lfor four representa- 
tive SEGUE stripes. As concluded in ^2.9.51 the Ry results for 
Ar < 1 are expected to be biased low. For A,- > 2, and outside 
the red giant region, the median Ry does not deviate appre- 
ciably from its canonical value. A more quantitative descrip- 
tion of this behavior is shown in Figure |33] For stars selected 
by 1 kpc < D < 2.5 kpc and Ar > 2.5 from / = 110° stripe, 
the median Ry is 2.90, with a mean of 2.95 and an rms of 
0.22 (determined from the interquartile range, the sample size 
is '--^9,000 stars). Given various systematic uncertainties that 
cannot be smaller than 0.1-0.2, as well as expected random 
errors ('-^O.l), the median Ry is consistent with the canoni- 
cal value of 3.1. We note that the width of the Ry histogram 
is about twice as large as the width of Ry determined using 
a fixed-Ry mock sample. Assuming that both widths are reli- 
able, which may not be strictly quantitatively true, the implied 
instrinsic scatter in Ry for / = 110° stripe is ^ 0.2. Results 
from other stripes are similar, with the median Ry showing a 
scatter of about 0.1. To illustrate this Ry variation. Figure |33] 
also shows the Ry distribution for I = 70° stripe, which has a 
median Ry of 2.80, and an rms of 0.15. We have tested for 
a possibility that the variation of the Ry distribution among 
stripes is due to calibration problems by comparing the me- 
dian residuals between observed and best-fit magnitudes for 
each SDSS camera column (twelve per stripe) and filter (in- 
cluding 2MASS filters). We did not find any evidence for 
photometric calibration errors larger than 0.01 mag. 



As shown in Figure |32] candidate red giant stars (top left 
corner) have consistently somewhat larger values of Ry (by 
about 0.2-0.4) than the typical star in the sample. Given 
that they are expected to be at much larger distances than 
main sequence stars, it is possible that they sample differ- 
ent types of dust. However, given fairly large range of lon- 
gitudes sampled by SEGUE stripes, this conclusion would 
imply that the dust in the Solar neighborhood (within 1-2 
kpc) has anomalously low Ry . A more plausible explanation 
for increased Ry is a bias due to slight differences in SEDs 
between red giants and main sequence stars. A preliminary 
analysis of the SDSS spectroscopic sample has revelead that 
spectroscopically-confirmed giants show an offset from the 
Covey et al. locus in the seven-dimensional SDSS-2MASS 
color space. Such an offset is, at least in principle, capable of 
inducing a bias in best-fit Ry. A detailed analysis of this bias 
and differences in SEDs between main sequence stars and red 
giants will be presented elsewhere. For the remainder of anal- 
ysis presented here, we simply exclude candidate red giant 
stars. 

A cross-section of the three-dimensional Ry map is shown 
in Figure [34] (recall that the Ry values are not reliable in re- 
gions of small Ar', see the rightmost panel for reference). For 
most of high-Ar regions, the median values are consistent with 
the canonical values. 

4. THE THREE-DIMENSIONAL DISTRIBUTIONS OF DUST AND 
STARS 

Best-fit stellar distance and extinction along the line of 
sight, Ar, determined here can be used to infer the three- 
dimensional distributions of dust and stars. The determina- 
tion of these distributions is not straightforward. In case of 
stars, complicated flux-color-extinction selection effects must 
be taken into account in order to obtain unbiased distributions. 
This analysis is best done with the aid of mock catalogs, such 
as those produced by Galfast. In case of dust, the complexity 
is further increased because the integral of dust volume den- 
sity along the line of sight is constrained, and not the density 
itself. To translate these constraints into a positive dust vol- 
ume density (more precisely, extinction per unit length as a 
function of position in the Galaxy), a careful statistical treat- 
ment of all errors and selection effects is mandatory. Since 
the full analysis is obviously far beyond this preliminary in- 
vestigation, we illustrate the potential of our dataset with two 
simplified analysis examples. 

4. 1 . The Spatial Distribution of Dust 

A coarse approximate map of the spatial distribution of dust 
in a given distance range can be obtained by subtracting two 
median A,- maps corresponding to the distance limits of the 
chosen range. This method is not statistically optimal, but it 
suffices for simple visualization. Figure |35] shows the result 
of such analysis for mean bin distances of 1.0, 1.5, 2.0 and 
2.5 kpc, with limiting bin distances 0.5 kpc larger and smaller 
than the mean distance. 

It is easily discernible that the dust structures observed at 
b ^2° and b ^ 13° are confined to 1-1.5 kpc distance range, 
while the structure seen at -3° < b <0° is due to dust at a 
distance of ^2.5 kpc and subtends < 1 kpc along the line 
of sight (an analogous panel for a mean distance of 3.0 kpc 
shows that this structure is mostly confined to smaller dis- 
tances). As discussed earlier, this ability to "bracket" dis- 
tances to dust clouds, and in turn to molecular clouds, is an 
important feature of our dataset. 
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Another projection of our dataset, the median as a func- 
tion of spatial coordinates, is shown in Figure [36] for all ten 
SEGUE stripes. Aside from the fact that data for each stripe 
also resolve the third direction (Galactic longitude), this pro- 
jection illustrates the integral constraint on the spatial distri- 
bution of dust. For each pixel, or a star in general case, the 
measured contains (noisy) information about the dust dis- 
tribution along the line connecting this pixel/star and the ob- 
serving point (the origin in this figure). With an appropri- 
ate model description of dust distribution, either parametric 
or non-parametric, these maps can be used to constrain 
the ni odel (for an example of similar analysis, see iJones et alJ 
1201 ll) . We point out that latitudes with most dust in a given 
stripe vary with the stripe longitude. For example, for the 
I ~ 130° stripe, the highest-extinction regions are found at 
positive latitudes, while for the / ~ 230° stripe, most dust is 
found at negative latitudes. 

4.2. The Spatial Distribution of Stars 

The spatial distribution of stars (the number volume den- 
sity) is shown in Figure|32] We have accounted for the change 
of volume with distance, but the variable distance limit due to 
faint flux cutoff and variable A,- is not taken into account and is 
clearly visible in the figure. To fully exploit these data for con- 
straining Galactic structure models, a three-dimensional dust 
map must first be derived from A,- constraints (or at least care- 
fully considered to mask high-extinction regions), and then 
one must apply color-dependent distance corrections. 

Nevertheless, several encouraging features are already dis- 
cernible in Figure|37] First, the sample seems fairly complete 
for distances below 1 kpc, corresponding to vertical distances 
from the plane of up to |Z| ^0.5 kpc. This volume is poorly 
explored by SDSS high-latitude data (e.g. see Figure 15 in 
JOS) and the dataset presented here will enable detailed stud- 
ies of the disk stellar number density profile for small |Z| (e.g., 
is the exponential profile valid within 100 pc from the disk 
mid-plane?). Second, the stellar number density at a fidu- 
cial location (say, at a distance of 0.5 kpc and Z = 0.3 kpc) 
significantly varies with Galactic longitude. This is expected 
behavior for an exponential disk profile in the galactocentric 
radial direction, and these data can be used to improve the ex- 
ponential scale length estimates for thin and thick disks (for 
more details, please see §4 in JOS). 

5. SUMMARY AND DISCUSSION 

This is the first analysis based on SDSS data that simultane- 
ously estimates intrinsic stellar color and dust extinction along 
the line of sight for several tens of millions of stars detected in 
the low Galactic latitude SEGUE survey. The fitting method 
and various assumptions are described in ^ Our main results 
are: 

1. The wavelength range spanned by the SDSS photomet- 
ric system and the delivered photometric accuracy are 
sufficient to constrain the intrinsic stellar SED and dust 
extinction along the line of sight. The minimum re- 
quired photometric accuracy of ^0.03 mag prevents 
non-unique solutions in most cases, and the accuracy of 
best-fit parameters scales roughly linearly with smaller 
errors. At the same time, this accuracy requirement ef- 
fectively limits the sample to about r < 19. 

2. Using the joint SDSS-2MASS photometry for stars at 
high Galactic latitudes, we confirmed the SDSS-based 



result from Sch2010 that the lO'Donne"lil(ll994l) redden- 
ing law can be rejected. We adopted the ICardelli et al.l 
(Il989h reddening l aw in this work, which is similar 
to the iFitzpatricM (Il999h reddening law adopted by 
Sch2010. Formally, both models are mildly inconsis- 
tent with the SDSS-2MASS data, but in practice pho- 
tometric implications of these differences are minor 
(^ 0.01 mag when A^ = 1). We recommend the coef- 
ficients listed in the first row in Table 1 for correcting 
SDSS and 2MASS photometry for interstellar dust ex- 
tinction. 

3. For stars detected by both SDSS and 2MASS, and when 
Rv is not a free fitting parameter, the best-fit intrin- 
sic stellar color and A^ for only-SDSS (four colors) 
and SDSS-2MASS (seven colors) fitting cases are sim- 
ilar Although SDSS samples reach much further than 
SDSS-2MASS samples at high Galactic latitudes (the 
distance limits for blue stars differ by about a factor of 
10), this is not the case at low Galactic latitudes because 
observed sources are much redder due to dust. The 
main benefit of only-SDSS samples is about a factor 
of 2 larger distance limit for blue main sequence stars; 
however, the limiting distance for red stars is smaller 
than for SDSS-2MASS case due to a necessary limit on 
the u band photometric errors. 

4. The SDSS photometry is not sufficient to reliably esti- 
mate Ry, with a realistic mock catalog implying errors 
of about 1. However, the addition of 2MASS photome- 
try significantly improves the accuracy of Ry estimates, 
with realistic mock catalogs implying errors of about 
0.3 when A,. ^ 1, and as small as 0.1 for A,. > 4. When 
Rv can be reliably estimated, we find that Ry = 3.1 can- 
not be ruled out in any of the ten SEGUE strips (at a 
systematics-limited precision level of ~ 0.1-0.2). Our 
best estimate for the intrinsic scatter of Ry in the re- 
gions probed by SEGUE stripes is < 0.2. 

5. Simultaneous fits for the intrinsic stellar SED and dust 
extinction along the line of sight allow for efficient 
recognition of candidate red giant stars in the disk. 
The selection method, which we dub "dusty parallax 
relation", utilizes the increase of dust extinction with 
distance, and identifies candidate giants as stars with 
anomalously large best-fit A^ for their best-fit main- 
sequence distance. 

6. The SDSS-2MASS photometi'ic dataset aflows robust 
mapping of the three-dimensional spatial distributions 
of main sequence stars and dust to a distance of about 
2 kpc (and Ar < 2 - 3). To extend this distance 
limits, deeper optical and infrared data are needed. 
With LSST and Wide-Field Infrared Survey Explorer 
(WISE) datasets (see below), the distance limit could 
be extended by close to a factor of 10. 

7. The three-dimensional spatial distributions of stars and 
dust can be readily analyzed with the datasets discussed 
here, which we make public (see Appendix B). For 
most scientific applications, we recommend tlie use 
of tlie SDSS-2MASS dataset witli fits based on aU 
seven colors, and with Ry fixed to its canonical value. 
For studies exploring the Ry variations, the use of the 
full SDSS-2MASS dataset is mandatory, with best-fit 
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Rv trustworthy only for stars with X^p^f < 2 and > 1 . 
Our fits represent a "stress test" for both SDSS and 
2MASS photometry, and we emphasize that careful 
quality cuts must be applied to avoid unreliable results! 



Given the results prese nted by Sch2010, iPeek & GravesI 
( 120101) . iJones et al] (1201 Ih and here, it is confirmed beyond 
doubt that there are some systematic problems with the nor- 
malization of SFD extinction map. Nevertheless, at high 
Galactic latitudes with small extinction these errors do not 
dominate over the photometric zeropoint calibration errors 
in SDSS data (0.01-0.02 mag), and at low Galactic latitudes 
most stars are embedded in dust and thus the SFD map is of 
limited use. 

Analysis described at the beginning of ^3.2| shows that the 
datasets analyzed here can robustly distinguish predictions 
made by the three popular models for the shape of dust extinc- 
tion curve. The O'Donnell model is clearly excluded, and the 
other two models do not provide a perfect fit to data either. On 
the other hand, the differences are very small and not much 
larger than systematic errors in photometry. The systematic 
photometric and other errors translate to a systematic uncer- 
tainty in Ry of about 0.1-0.2. We did not detect any devia- 
tions from the canonical value = 3.1 at this precision level. 
We reach the same conclusion bv Jones et al. (20 11), but here 
we obtained several times smaller errors due to a much wider 
wavelength range of utilized photometry. This uniformity of 
dust properties within a fairly large volume (distance limit of 
the order 1 kpc) probably implies that the ISM dust is well 
mixed during its lifetime (Draine 201 1). 

Last but not least, it will be very informative to di- 
rectly compare the results presented here with those ob- 
tained by ot her methods, such as near-infrared color ex- 
cess method (Lomb ard! & Alvelll2001l: [Lombardi et al.' 201 It 
[Majewski et al. 2011]), Ha-based method (Saleetal. 20091 
and Wolf method that is sensitive to gray dust ( Yasud a et al] 
12007; Gorbikov & Brosch.201 0). to uncover and quantify var- 
ious systematic errors that are likely to exist in all methods. 

The results presented here will be greatly improved by 
several upcoming large-scale , deep optical s urveys, includ- 
ing the Dark Energy Survey dFlaughed 120081). Pan-S TARRS 
dKaiser et al.l |2002), and the LSST dlvezic et al.1 l2008bl) . 
These surveys will significantly extend the faint limit of the 
sample analyzed here (in case of LSST by 5 mag) and 
are likely to provide more reliable photometry due to multi- 
ple observations and the use of photometric methods designed 
for crowded fields. Although 2MASS is too shallow to fully 
complement these new optical surveys, it will still provide 
very useful constraints in high-extinction region s. Further- 
more, the recently released WISE data (Wright et al.ll2010h 
will provide supplemental constraints with its Wl band at 3.4 
/im, which reaches about 2 mag deeper than 2MASS K band. 



These new datasets are thus certain to provide valuable new 
information about the dust and stellar distribution within the 
Galactic disk beyond the current limiting distance of a few 
kpc. 
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APPENDIX 
SQL QUERY EXAMPLE 

The following SQL query was used to select and download data for all SDSS stars with spectroscopic and proper-motion 
measurements (see http://casjobs.sdss.org/CasJobs). 

SELECT 

round (p . ra, 6 ) as ra, round (p . dec, 6) as dec, 

p. run, p.camcol, p. field, comments are preceded by 

round (p . extinction_r, 3 ) as rExtSFD, r band extinction from SFD 

round (p . modelMag_u, 3 ) as uRaw, N.B. ISM-uncorrected model mags 

round (p . modelMag_g, 3 ) as gRaw, rounding up 
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round (p . modelMag_r , 3 ) as rRaw, 
round (p . modelMag_i , 3 ) as iRaw, 
round (p . modelMag_z , 3 ) as zRaw, 
round (p . modelMagErr_u, 3 ) as uErr, 
round (p . modelMagErr_g, 3 ) as gErr, 
round (p . modelMagErr_r , 3 ) as rErr, 
round (p . modelMagErr_i , 3 ) as iErr, 
round (p . modelMagErr_z , 3 ) as zErr, 

(case when (p. flags & '16') = then 1 else end) as ISOLATED, 

ISNULL (round (t .pmL, 3) , -9999) as pmL, proper motion data are set to 

ISNULL (round (t .pmB, 3) , -9999) as pmB, -9999 if non-existent (NULL) 

ISNULL (round (t .pmRaErr, 3) , -9999) as pmErr if pmErr < no pm data 

INTO mydb. dustSample 

FROM phototag p LEFT OUTER JOIN propermotions t ON 

(p.objID = t.objID and t. match = 1 and t.sigra < 350 and t.sigdec < 350) 
quality cut on pm 

WHERE 

p. type = 6 and select unresolved sources 

(p. flags & '4295229440') = and '4295229440' is code for no 

DEBLENDED_AS_MOVING or SATURATED objects 

p. mode = 1 PRIMARY objects only, which implies 

IBRIGHT && (IBLENDED | | NODEBLEND | | nchild == 0)] 

p.modelMag_r < 21 adopted faint limit 

the end of query 

DATA DISTRIBUTION 

All data files, as well as a detailed description of their content, are available from a public data repositor}0 Due to the large 
data volume, we separate our catalogs into four groups. We fit stellar SEDs twice for all 10 SEGUE strips: once with selective 
extinction fixed at Ry = 3.0, and a second time with Ry as a free fitting parameter (limited to the range 1 -7.9). Similarly, we 
present the only-SDSS and SDSS-2MASS datasets separately. For the Ry = 3.0 case, the data files in each dataset (only-SDSS 
and SDSS-2MASS) are defined by Galactic coordinates, and are designed to contain fewer than 10 million stars each. For the 
fiee-Ry case, we distribute only the data from SEGUE strips with \b\ < 30° because Ry is poorly constrained at higher Galactic 
latitudes with small extinction. This data organization allows users to download data for a relatively small region of sky without 
the burden of downloading the whole dataset. These datasets are made available in two formats: as FITS tables, and as plain 
ASCII text files. 

All of the data files contain SDSS astrometry and photometry (and proper motions), the SFD value for Ar, and best-fit model 
parameters (including a best-fit distance estimate). Additionally, the SDSS-2MASS data files also contain 2MASS astrometry 
and photometry, the only-SDSS best-fit parameters, and the SDSS-2MASS best-fit parameters. 

We emphasize that our fits represent a "stress test" for both SDSS and 2MASS photometry, and thus careful quality cuts must 
be applied to avoid unreliable results! 

DISCUSSION OF THE METHODOLOGY 
Here we provide a more detailed discussion of two aspects of methods discussed in ^ 



Closing the System of Equations 

Stellar colors constrain reddening due to dust, e.g., Aug = A„-A^, rather than dust extinction, here A„ and Ag. Therefore, 
when inferring the amount of dust extinction, both in case of single stars that are projected onto the unreddened stellar locus 
in the multi-dimensional color space, and in case of color offsets of the whole stellar locus at high Galactic latitudes, there is 
always one constraint /ewer than the number of photometric bands. A convenient way of thinking about this "missing" equation 
is that dust extinction is described by its "scale" Ar and four (or seven in SDSS-2MASS case) measures of the scaleless shape 
of the extinction curve, C\ = AxjA^. Three different approaches can be used to "close" this system of equations, and to break 
"reddening-extinction" degeneracy (we do not discuss the best approach, based on known distance modulus, DM, and absolute 
magnitude, M, , which directly constrains Ar via r = Mr + DM+Ar, because our dataset does not include a parallax distance to the 
vast majority of stars). 

The first approach assumes that A,, is provided as an additional input, for example, from the SFD map as A^^^. In this case, 
Aug = (Cu - Q) Af™, and it is easy to show that 

C«=l + ^5^ (CD 

http://www.astro.washington.edu/users/ivezic/r_datadepot.html 
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C. = l-^- (C4) 

If there are systematic errors in Af'^'^, they will be propagated to C,„. Such effects can be tested for by tracing the variation of 
resulting Cm across the sky, and by correlating deviations with A^^^. In particular, given many lines of sight, it is possible to fit a 
spatially invariant model for errors in A^'^^ (e.g., an additive and a multiplicative error). 

To illustrate the impact of errors in A^^^ on C,„ determined with this method, we computed Acoior using the CCM model with 
Rv = 3.1 and true Ar = 1, and we assumed a multiplicative error in Af'^'^. A correction factor of 0.95 produces the overall best- 
fit Ry = 3.28, with the / and z band constraints on Ry biased to even higher values. For a correction factor of 0.9, the best-fit 
Ry = 3.48, with the / and z extinction values barely consistent with the CCM extinction curve. For additive errors such that the 
true Ar = A^™ + 0.05, the best-fit Ry varies from 2.20 forA^'"" = 0.1, to 2.93 at for A^™ = 1. Therefore, this method is quite 
sensitive to systematic errors in Af'^'^ and should be used with caution. 

The second approach uses a model -based extinction curve to predict Cm as a function of single parameter Ry. Given that there 
are three "spare" constraints, model predictions can be tested for self-consistency (and perhaps used to select the "best" model). 
This method results in estimates for A,, and can be used to test external maps, such as SFD, though only in a model-dependent 
way. 

The third approach, adopted here, is to assume (fix) one value of C,„ and solve for A^ and all remaining C,„. While at first this 
approach sounds arbitrary, it becomes sound when SDSS data are augmented with 2MASS data. The reason is that the effective 
wavelength for 2MASS K band is 2.2 /xm, which in this context is almost as large as infinitjQ When using both SDSS and 
2MASS, A,- is estimated using offsets of the r-K color distribution. The main reason why this approach works is the fact that 
AK/Ar is small (0.132) and varies little with Ry, and among all plausible dust extinction models. For example, for Ry > 2 all 
models predict variations of AxjAf not exceeding 20%. This variation translates to only about a 3% error when the r-K shift 
is interpreted as A,--Ak = (1 -Ajf/A,.)A,- = 0.868Ar. When using this approach, there are seven colors constructed with eight 
photometric bands, and the result is estimates for Ar and six Cm. 



Methods for Quantifying Color Offsets for the Stellar Locus 

What is the optimal method for measuring A„g, Ag^, etc., using the stellar locus? If we think of the stellar locus in a two- 
dimensional (2D) color-color diagram as of an "image", then we essentially "slide" the image of the reddened sample to perfectly 
align with the image of the "intrinsic" dereddened locus. This alignment can be performed in each 2D color-color diagram, or 
alternatively all four color shifts can be determined simultaneously in the 4D color space. At the other extreme, the color shifts 
can be determined using ID projections of each color, as in the "blue tip" method proposed bv lSchlaflv et al.l (12010'). 

If there were no astrophysical systematics and measurement error distributions were fully understood, these methods should 
produce identical results (e.g., the g-r offsets estimated from the g-ry&. u-g and r- i vs. g-r diagrams would be statistically 
consistent). However, there are astrophysical systematics, such as distance, age and metallicity effects, that may introduce various 
biases. For example, M dwarf stars in SDSS sample can be as close as 100 pc and thus be within dust layer, and the "blue tip" 
is sensitive to age and metallicity of turn-off stars that define it. The idea behind the principal colors method employed here is 
to avoid distance effects by considering only stars bluer than M dwarfs, and to mitigate age and metallicity effects by measuring 
shifts perpendicular to the locus. The reason for the latter is that age variation "extends" or "shortens" the locus (i.e., shifts the 
"blue tip"), but does not strongly affect its position in the perpendicular (Pj) direction. When considering metallicity, systematic 
effects are a little bit more complicated, but are mostly confined to the u band. For blue star s, the g-r color is essentially a 
measure of effective temperature with negligible dependence on metallicity jlvezic et al.ll2008ah . At a g iven g-r color, the u — g 
color depends on metallicity (it becomes bluer as metallicity decreases, see the top right panel in fig. 2 in llvezic e t al. 2008a). For 
example, at g-r=0.3, the u-g color varies by about 0.2 mag as the metallicity varies from the median thick-disk value (-0.5) to 
the median halo value (-1.5). This shift is not parallel to the locus in the g-r vs. u-g color-color diagram, so it does have some 
effect on the P2 distribution. However, already at g-r = 0.5, the fraction of halo stars in SDSS sample is sufficiently small that this 
effect becomes negligible (because such red halo stars are too faint to be detected by SDSS). Hence, in the range 0.5 < g-r < 1 .2, 
only the dust reddening (and photometric calibration errors, of course!) can significantly shift the locus perpendicularly to its 
blue part (even in the u-g vs. g-r diagram). An added benefit from the signal-to-noise ratio viewpoint is that the P2 distributions 
are very narrow, an advantage that mitigates the fact that the reddening vectors are measured only along P2 directions. 
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TABLE 1 

Observational Constraints and Model Values for the Extinction Curve, Cx = Ax/A,- 



Region 


M 


g 






./ 


H 


K, 


S82 


1.810 


1.400 


0.759 


0.561 


0.317 


0.200 


0.132 


North 


1.750 


1.389 


0.750 


0.537 


0.297 


0.180 


0.132 


CCM 


1.814 


1.394 


0.764 


0.552 


0.327 


0.205 


0.132 


F99 


1.795 


1.415 


0.748 


0.554 


0.308 


0.194 


0.132 


OD 


1.813 


1.406 


0.783 


0.562 


0.325 


0.205 


0.132 



Note. — The first two rows list observational constraints for the shape of the extinction curve, Cx = Ax/ A,-. The value of Cx in the K band was assumed to 
he 0.132. The first row corresponds to the so-called SDSS Stripe 82 region (defined by 300° < R.A. < 60° and \Dec\ < 1.27°), and the second row to a northern 
region defined by 30° < h < 45° and 0° < / < 10°. The last three rows list model predictions computed for an F star spectrum and the best-fit value oi Ry 
(CCM=Cardelh et al. 1989: Rv = 3.01; F99=Fitzpatrick 1999: Rv = 3.30; OD=0'Donnell 1994: Rv = 3.05). 



TABLE 2 

Adopted Extinction Coefficients, CaCRv) 



Rv 




u 


g 


(' 


z 


J 


H 


K, 


Source 


2.0 


2 


280 


1.579 


0.702 


0.453 


0.264 


0.166 


0.107 


CCM 


2.5 


1 


998 


1.467 


0.740 


0.513 


0.302 


0.190 


0.122 


CCM 


3.0 


1 


817 


1.395 


0.764 


0.552 


0.326 


0.205 


0.132 


CCM 


3.1 


1 


788 


1.384 


0.768 


0.558 


0.330 


0.208 


0.134 


CCM 


3.1 


1 


855 


1.446 


0.743 


0.553 








Sch2010 


3.1 


1 


857 


1.439 


0.725 


0.517 


0.250 


0.131 


0.068 


F99 


4.0 


1 


598 


1.308 


0.793 


0.598 


0.356 


0.224 


0.144 


CCM 


5.0 


1 


470 


1.257 


0.810 


0.625 


0.373 


0.234 


0.151 


CCM 



Note. — An illustration of the dependence of the adopted extinction curve, Cx = Ax/Ar on Rv (Cr = 1 by definition; see also Figure[8). The second line with 
Rv = 3.1 lists the values suggested by Schlaflv & Finkbeiner 2010, and the third fine with Rv = 3.1 lists the values computed using eq. 5 from Fitzoatrick & Mass3 
12009 with a = 2.50 (constrained by C,- = 1), and using Av/A^ = 1.20 and the effective wavelengths from tSchlaflv & Finkbeine r 2010) for the SDSS bands, and 
1.25 fim, 1.65 fim, and 2. 17 /im for the 2MASS JHK bands, respectively. Both fines are presented for a comparison with the adopted CCM extinction curve. 
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r-K 



Fig. 2. — The distribution of unresolved SDSS sources witli 2MASS detections in tfie A — r vs. r — K color-color diagrams, witli A = u,g.i,z,J and H. Tlie 
source density is sliown as color-coded maps, and it increases from black to green to red. The two arrows marked PCj and PCt in the top left panel illustrate the 
"principal color" axes discussed in text and u sed to track the locu s shifts due to interstellar dust reddening. The dashed vector in each panel shows the reddening 
vector for A, = 2 and standard Rv = 3.i dust <CardeUi et al.ll989l) . 



22 



Berry et al. 




Fig. 3. — Model predictions for tiie extinction curve siiape as a function of for tliree different models: O'D fO'Donnell"1994'), F99 fFi tzpatrick|[T999l) . and 
CCM (Cardelli et al. 1989), evaluated for stars with three different effective temperatures (as listed in the legend, in Kelvin). The left panel shows Cx =Ax/Ar 
for A = r, i,z,J,H,K) (top to bottom, respectively); the right panel is analogous, except that the ratios based on colors (u-g, g-r,r-i,i-z, z—J,J—H, and 
H—K) are shown. As expected, most of the sensitivity to Rv comes from the blue bands (u and g). 




Fig. 4. — A comparison of the constraints on the extinction curve shape (the three plus symbols, with approximate Icr uncertainty limits shown as ellipses) and 
three model predictions (see Figure[3]for legend; the three crosses along the curves correspond to Rv=2.6, 3.1, and 3.6). The pink symbol corresponds to the 
Stripe 82 region (southern Galactic hemisphere), the brown sy mbol to the northern Galactic hemisphere, and the blue symbol is the constraint from t he Schlafly 
et al. (2010) analysis. The blue (Fitzpatrick 1999) and green (Cardelli et al. 1989) models are in fair agreeement with the data, while the red model 1(3 ' Donnell 
119941) predicts unacceptable values of the (Ar— A,)/(A,— A,) ratio for all values oi Ry (see also Figure|5). 
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Fig. 5. — Constraints on Ry based on the CCM jCardelli et al.ll 19891) dust reddening law. Only the SDSS bands, which provide the strongest constraints on 
Ry. are shown (see the legend). The dashed line shows the overall constraint on Ry (posterior probability distribution for a flat prior), with the best-fit value of 
=3.01 ±0.05. 
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Fig. 6. — Analogous to Figure|5] except that lO'Donnell fl"994l) dust reddening law is used. The predicted extinction in the i band is inconsistent with constraints 
from other bands. 
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Fig. 7. — The left panel shows the CCM (black, solid line) and the O'Donnell (red, dashed Hne) dust reddening laws as function of wavelength. The right panel 
displays the fractional difference between the two dust extinction models with the SDSS filter transmission regions overlaid (vertical, blue, dashed lines). Given 
the filter transmission regions, the largest integrated difference is expected in the i band, which is what we observe. 




Fig. 8. — The adopted Ax/ A,- ratio, shown as a function of Rv. for A = (ugrizJHK), from top to bottom (Ar = 1). The curves are computed for an F star using 
the CCM t'Carde lli et al. 1989) dust reddening law. 
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Fig. 9. — An illustration of the constraints on intrinsic stellai' colors, extinction in the r band, Ar, and the ratio of total to selective extinction, R\/ . In both 
diagr ams, the line arly-spaced contours show the main stellar locus as observed at high Galactic latitudes. The dashed lines mark the median stellar locus 
from ICovev et al J j2007). In the left panel, the dot marked "Obs" represents a hypothetical observation. Depending on the adopted Ry, as marked, different 
combinations of intrinsic stellar colors (i.e., the position along the stellar locus) and A, are consistent with the observed g-r and r—i colors. Multiple solutions 
are possible even for a fixed value of R\i. The three solutions marked 1-3 con'espond to (Rv,A,)= 1:(2.2,1.0), 2:(5.0,2.2), and 3:(5. 0,6.0). As shown in the right 
panel, these degeneracies are broken if the i — color is also available: the three {RyJ^r) combinations have different reddened i — colors which breaks the 
degeneracy between the intrinsic stellar color and A, . The degeneracy is broken because the reddening vectors in the right panel are nearly parallel despite very 
different Ry values. 




Fig. 10. — A Monte Carlo study of best-fit stellar model errors (parametrized by the g-i color) as a function of photometric errors, for a fiducial star with 
g — i=l.95 andAr=1.5 (the abscissa, Ag-(=true - fit). The photometric errors are generated from Gaussian distributions with widths equal to 0.01 mag (top left), 
0.02 mag (top right), 0.04 mag (bottom left) and 0.08 mag (bottom right). The errors in the best-fit g-i aie about twice as large as assumed photometric errors. 
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Fig. 1 1 . — Analogous to Fig. 1101 except that the errors in the best-fit Ar are shown (AA,)=true - fit). Note that for large photometric errors (the bottom two 
panels), the A, error distribution becomes bimodal; the additional mode corresponds to a solution with a bluer star with more reddening. 
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Fig. 12. — Analysis of the covariance in the best-fit values for A,- and g — i using a simulated dataset. The panels show the distributions of the best-fit values for 
Ar and g-i for two different fiducial stars (left column: a blue star with true g-i=OA; right column: a red star with true g — (=3.0), and two different extinction 
values (top panels: A,- = 1; bottom panels: Ar=3). Photometric eiTors in the ugriz bands are generated using Gaussian distributions with (T=0.02 mag (uncorrelated 
between different bands). Note that the A,- vs. g-i covariance is larger for the blue star, and does not strongly depend on assumed A,-. 



Milky Way Tomography IV: Dissecting Dust 3 1 




Fig. 13. — Analysis of a Galfast simulated SDSS-2MASS sample from a SEGUE strip (/ ~ 110° and \b\ < 5°). The r vs. g-i color-magnitude diagrams in 
the top two panels explain why the fraction of giants is much larger than observed at high Galactic latitudes. The same simulated sample, defined by obsenvd 
extincted magnitude cuts 14 < r < 21 and K < 14.3 (Vega) is shown in both panels. The left panel is constructed using im-extmcted magnitudes, and the right 
panel with "observed" magnitudes (note the offset of the y axis by 2 mag). The horizontal dashed line in the left panel shows the SDSS saturation limit; stars 
above this line are dominated by red giants (the "plume" towards towards g-i ~ 1). The diagonal dashed fine shows the magnitude limit for main sequence stars 
with K < 14.3 and no dust extinction. The reddening arrow corresponds to Ar = 2 and Rv = 3. 1 CCM extinction curve. The bottom left panel is analogous to the 
top left panel, except that the SDSS-based best-fit values for g—i and A,- are used. The bottom light panel shows the difference between the input value of Ar and 
the best-fit values, as a function of the latter. The dashed line is added to guide the eye. The root-mean-square scatter for the A, difference (rms for y axis) is 0.33 
mag, and the bias for large A, is about 15% (an overestimate of A,, due to color-A,- degeneracy, see text). 
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Fig. 14. — A comparison of six SDSS-2MASS color-color diagrams for data from the SEGUE / ~ 110° strip (the top six panels), and for a mock catalog 
produced with the Galfast code (the bottom six panels). The color-coded contours shows the source counts on a linear scale. The two dashed arrows show 
reddening vectors for A,- = 2 and Ry = 2 and 4. The locus of circles shows the Covey et al. 1 2007) empirical SED library and illustrates the morphology of the 
same diagrams observed at high Galactic latitudes (and corrected using the SED map; typically A,- ~ 0.1). The two sets of diagrams are encouragingly similar, 
with a few detailed differences: the observed diagrams have more outliers, and a few diagrams (e.g., J-K vs. i — z and i-z vs. r-i) imply different reddening 
vectors than used in simulations (Ry = 3.1). 
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Fig. 15. — A comparison of the best-fit Ry values for SDSS-2MASS (narrow histogram) and only-SDSS (broad histogram) cases, using a simulated Galfast 
mock catalog. The input value is fixed to i?v =3.1. The equivalent Gaussian widths determined from the interquartile range are 0.1 and 1.2, respectively. 
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Fig. 16. — The left panel shows the relationship between best-fit Ar and distance computed usin g th e S PSS-based best-fit stellar color and a photometric 
parallax relation appropriate for main sequence stars, for the same simulated sample as in Figures [T3l and [T4l The color-coded map shows the counts of stars on 
a Hnear scale. The dashed lines isolate candidate red giant stars that have small distances and large A,-. The top right panel is analogous to the top left panel, 
except that the best-fit values correspond to SDSS-2MASS data. The bottom left panel shows the median input absolute magnitude (M,-) for stars in each pixel, 
color coded according to the legend in the lower right corner. The "red giant region" in the top two panels is dominated by gia nts ( M,- < 3). The bottom right 
panel shows the difference between true absolute magnitude and an estimate obtained from the "dusty parallax relation" (see eq. |13K as a function of the latter. 
The dashed line has is added to guide the eye. The root-mean-square scatter between the two magnitudes (rms for y axis) is 1 . 1 mag. 
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Fig. 17. — The distribution of the best-fit Xpjf (Rv = 3.1), with differential distributions in the left two panels, and cumulative distributions in the right two 

panels. The top two panels con'espond to the north Galactic cap region (h > 45°) and the bottom two panels to the SEGUE / ~ 1 10° strip, limited to \b\ < 5° (a 
high-extinction region). The solid lines are used for SDSS-only fits, and the dashed lines for fits to SDSS-2MASS data. In the two left panels, the top solid line 
corresponds to subsamples of stars with r < 20, and the bottom solid line to stars with 20 < r < 21. The top dashed line corresponds to the full SDSS-2MASS 
sample, and the bottom dashed line to subsamples with K < 13.9 (Vega scale, approximately con'esponding to K band eiTors up to 0.05 inag). The solid lines 
in the right panels correspond to the full SDSS sample, and the short-dashed lines to the full SDSS-2MASS sample. The dot-dashed and long-dashed lines 
correspond to Xpjf distributions with 2 and 5 degrees of freedom. The long-dashed line in the top left panel corresponds to the Xpjf distribution with 5 degrees of 
freedom for the full SDSS-2MASS sample. 
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Fig . 18 . — Analysis of the best-fit results for A,- in the low-extinction region with h > 30° . The top left panel shows the median Ar in 0.6 deg^ pixels in Lambert 
projection. The values are linearly color-coded according to the legend. Stars with 15 < r < 18 and < 2 from only-SDSS sample with fixed Ry are used for 
the plot. The median difference between Ar and the values given by the SFD map are shown in the top right panel. Note the striping reminiscent of the SDSS 
scanning pattern. The bottom left panel shows the median difference between observed and best-fit model magnitudes in the band. The bottom light panel is 
analogous to the top right panel, except that only the subset of stars also detected by 2MASS (K < 14) a nd with full SDSS-2MASS fits (fixed Ry) are used. Note 
the much better agreement with the SFD values than in the top right panel. For more details, please see ^3.1.21 
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Fig. 19. — The color-coded maps show the best-fit A,- based on SDSS data for the ten analyzed SEGUE stripes. Each stripe is limited to the range of |i>| < 5°. 
A fixed Rv = 3.1 is assumed. The legend above each panel shows the color scale, and each 12x 12 arcmin- pixel shows the median Ar. For each stripe, three 
distance ranges are shown: 0.3-0.6 kpc (left), 1-1.5 kpc (middle) and 2-2.5 kpc (right). It is assumed that all stars are on main sequence when estimating distances. 
Only stars with best-fit Xpjf < 2 and outside the red giant region (selected here by Ar < 1.5 + 1.5Dkpc) Bie used for the plot. The top right panel shows the sky 
coverage of the full analyzed dataset. 
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Fig. 20. — The counts of stars in the only-SDSS case best-fit Ar vs. best-fit main sequence distance diagram for four SEGUE strips (top left: / = 70°, top right: 
/ = 110°, bottom left: / = 150°, bottom right: / = 187°; for all panels \b\ < 5°). Only stars with Xpdf < 2, r < 19 and error in the u band below 0.05 mag are 
used. Counts are normalized to the maximum value and color coded on the same linear scale, from blue (low) to red (high). The two dashed lines mark a region 
dominated by red giant stars (the top left corner). 
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Fig. 21. — Analogous to Figure llOl except for best-fits based on SDSS-2MASS sample (only stars with Xpif < 2 and < 15 on Vega scale are used). Note the 
larger fraction of red giant stars in the top left comer, and a smaller distance limit, compared to Figure |20l and that the fraction of giants decreases with Galactic 
longitude. 
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Fig. 22. — A comparison of six SDSS-2MASS color-color diagrams using data from the SEGUE / ~ 1 10° strip (the top six panels; same as the top six panels 
in Figure [14] except that here only st ars with Xpjf < 2 are used), and the best-fit model colors based on SDSS-2MASS dataset (the bottom six panels, in the 
same order). The thick lines show the lCovev et alj (2007) empirical SED library and illustrate the morphology of the same diagrams observed at high Galactic 
latitudes. The two sets of diagrams are encouragingly similar: fits to intrinsic stellar SED and dust extinction on per star basis are capable of reproducing the 
morphology of observed diagrams in highly dust-extincted regions. 
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Fig. 23. — A compaiison of the best-fit g — i (left panel) and Ar (right panel) values obtained with a fixed = 3. 1 for SDSS-2MASS sample from the SEGUE 
/ ~ 110° stripe, using two different fitting methods. The abscissae show the best-fit values obtained using only-SDSS dataset (four fitted colors), and the ordinates 
correspond to the residuals of the SDSS-2MASS (seven fitted colors) minus the only-SDSS datasets. The number density of stars increases linearly from black 
to blue to red. The dashed lines are added to guide the eye. 
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Fig. 24. — Analysis of the differences in the best-fit A,- between fits based on SDSS-2MASS dataset (first panel from the left) and those based on only-SDSS 
data (second panel). Only stars with best-fit X^pjf < 2, r < 20 and main-sequence distance 0.5 — 1 kpc are used for the plot. The top legend shows the coloring 
code for these two panels, and each 6x6 arcmin^ pixel shows the median Ar for stars with Xp^f < 2. The third panel shows the median difference between the 
two best-fit A, values (the second panel minus the first panel), with the color coding using the same palette, except that the limits are itO.l mag. The fourth panel 
shows the median difference in A, (i.e., the third panel) for 0.2° wide bins of Galactic latitude. 
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Fig. 25. — The color-coded maps show the best-fit A,- based on only-SDSS dataset for the SEGUE / ~ 110° strip. The legend shows the color scale, and each 
12x12 arcmin^ pixel shows the median A,-. The three panels correspond to main-sequence distance range: 0.3 — 0.6 kpc (left), 1 — 1.5 kpc (middle) and 2 — 2.5 
kpc (right). Only stars with best-fit X^,df ^ ^ outside the red giant region (selected here by < 1 .5 4- 1 -SDy^p^ ) are used for the plot. 




Fig. 26. — Analogous to Figure l25] except using the SDSS-2MASS dataset and different distance slices (left to right: 0.1—0.5 kpc, 0.5-0.7 kpc, 0.7-0.9 kpc, 
0.9- 1.1 kpc). Note the abrupt increase in A, for stars towards h ^ +2° that are more distant than 0.9 kpc. 
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Fig. 27. — Analysis of the differences between best-fit A,- values (left panel, based on only-SDSS data; SDSS-2MASS version looks similar) and the SFD 
values (second panel) for stars with X^,df < 2 and main- sequence distance in the range 0.8 — 1.2 kpc. The third panel shows the difference of the two A, values 

(the second panel minus the first panel). Each 6x6 arcmin- pixel in the first three panels is color coded according to the top legend. The fourth panel shows the 
median best-fit Ar (blue line) and the median SFD value (red line) for 0.2° wide bins of Galactic latitude. If the SFD maps are correct, then the dust structures 
discernib le in the two right panels at fo ~ 0° and h ~ +2° must be more distant than ~1 kpc. This conclusion is independently confirmed for the latter dust cloud 
in Figure l26l 
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Fig. 28. — The bottom four panels show the \h\ < 5° subregion of the panels shown in Fig. 1271 The top three panels show the mid-IR (left), CO (middle) 
and radio continuum (right) maps on approximately the same scale (obtained using "The Milky Way Explorer" by Kevin Jardine). The few small irregular 
white regions in the bottom three maps do not contain any stars with good photometry. Assuming that the SFD map is not grossly incorrect, the dust extinction 
determined here implies that most of the molecular cloud structures seen in the top middle panel must be more distant than ~1 kpc. 
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Fig. 29. — A comparison of the best-fit g — i (left panel) and A,- (right panel) values for two different treatments of Ry, for stars in the 1=1 10° SEGUE strip 
(using SDSS-2MASS data). Only stars with best-fit X^,jf < 2, r < 20 and K < 13.9 (Vega) are used for the plot. The abscissae show the best-fit values obtained 
for a fixed Ry = 3.1 and the ordinates correspond to the residuals of the differences in the best-fit values when Ry is treated as a free-fitting parameter ("free Ry" 
- "i?v=3.1"). The number density of stars increases from black to blue to red. The dashed lines are added to guide the eye. 
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Fig. 30. — The Ar map analogous to Figure 1211 except that Rv is treated as a free pai'ameter. Only stars with \h\ < 5°, Xpdf < 2, r < 19, and K < 15 (on Vega 
scale) are used. 
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Fig. 3 1 . — A comparison of three different types of best-fit SEDs: using only SDSS data witli fixed Ry = 3.1 (blue line), and using joint SDSS-2MASS dataset 
with fixed Ry (green line) and with free Ry (red line). As demonstrated by the similarity of best-fit lines, the differences in best-fit parameters, listed in each 
panel, are due to degeneracies between intrinsic stellar color, amount of dust and Ry. The shown cases correspond to blue and red stars (top row vs. bottom row), 
and small and large Ar (left column vs. right column). 
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Fig. 32. — Similar to Figure 1211 except tliat Ry is treated as a free parameter, and tlie color-coded map shows the median value of Ry (ranging from blue for 



= 1.5 to red for Ry = 4.5, green corresponds to Ry = 3; see the legend in the bottom right panel). Only stars with \b\ < 5° 



\.pdf 



< 2, r < 19, and K < 15 



(on Vega scale) are used. Note that red giant stars (top left corner) have consistently larger values of Ry, and that consistently Ry < 3 when A, < 1 for main 
sequence stars. In other regions in this diagram where Ry is determined robustly, Sv = 3.1 cannot be ruled out in any of the ten SEGUE stripes at a precision 
level of ~ 0.1 -0.2. 
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Fig . 33 . — A comparison of the best-fit Ry values for SDSS-2MASS free-Sv case and stars with distances in the 1 .0-2.5 kpc range and A,- > 2.5, selected from 
/ = 70° (red, left histogram) and 1 = 1 10° (blue, light histogram) stripes (other selection criteria are the same as for stars plotted in Figure[32). 
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Fig. 34. — The first three panels show the median Ry obtained using SDSS-2MASS sample for the SEGUE / = 1 10° strip, and for distance range 0.5—0.7 kpc 
(first), 0.7-0.9 kpc (second), and 0.9-1.1 kpc (third). Only stars outside the "red giant" region, see FigureOTI and with X^,df < 2, r < 21 and K < 14.3 (Vega) 

are used for the plot. The pixel size is 6x6 arcmin^, and the Ry coloring scheme is shown at the top. The fourth panel shows for reference the best-fit A, , for the 
distance slice 0.9- 1.1 kpc. Note that Rv is not reliable for A,. < 2 (black and blue regions in the fourth panel). 
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Fig. 35. — Illustration of the three-dimensional dust distribution for SEGUE stripe / ~ 110° at mean distances of 1.0, 1.5, 2.0 and 2.5 kpc, using only-SDSS 
sample and dxed-Ry fits. Unlike other figures that show the median Ar along the line of sight, this figure shows the differences in the median A,- (per 12x 12 
arcmin^ pixel) for samples at distances between the quoted distance and limiting distances 0.5 kpc larger and smaller than the mean distance (e.g., the first panel 
shows the difference between the median A,, for 0.5-1.0 kpc and 1.0-1.5 kpc subsamples). It is easily discernible that the dust structures observed at ~ 2° and 

~ 13° are confined to 1-1.5 kpc distance range, while the structure seen at -3° < i < 0° is due to dust at a distance of ~ 2.5 kpc (an analogous panel for a 
mean distance of 3.0 kpc shows that this structure is mostly confined to smaller distances). Note that the Hnear extent perpendicular to the line of sight of a given 
angular size is 2.5 times larger in the last than in the first panel. 
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Fig. 36. — The median best-fit A, (extinction along the line of sight) is shown as a function of distance from the Galactic plane, Z, and distance along the 
plane, D„, for 10 SEGUE stripes (this is not a cross-section of three-dimensional dust distribution!). The best-fit Ar are based on the SDSS-2MASS dataset and 
fixed-Rv fitting case, for stars with XJ,jf < 2 and K < 15 (Vega). Each pixel is 50x50 pc^ and subtends 2.5 deg wide stripe in the perpendicular (longitude) 
direction. The color scheme increases linearly from blue to red with a varying maximum value; 5 for the first three panels, 4 for the next three, and 3 for the last 
four panels. 
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Fig. 37. — The local volume n umbe r density of stars is shown as a function of distance from the Galactic plane, Z. and distance along the plane, Dry. for the 
same samples as shown in Figure [36l The color scheme shows the counts on log scale with the same arbitrary normalization for all stripes. The fall-off of the 
stellar volume number density at distances beyond ~ 1 kpc is due to the stellar color-dependent sample distance limit and does not reflect the disk structure. Note 
the variation of counts with Galactic longitude (the top four panels are closer to the Galactic center and contain more stars per unit volume). 



